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Ph. Abstract: Recently, Berenstein et al. have proposed a duality between a sector of 

M = 4 super- Yang-Mills theory with large R-charge J, and string theory in a pp-wave 
background. In the limit considered, the effective 't Hooft coupling has been argued 
to be A' = gY M N/J 2 = l/(fip + a') 2 . We study Yang-Mills theory at small A' (large 
/i) with a view to reproducing string interactions. We demonstrate that the effective 
genus counting parameter of the Yang-Mills theory is g\ = J 4 /N 2 = (47rg s ) 2 (^p + a f ) 4: , 
the effective two-dimensional Newton constant for strings propagating on the pp-wave 
background. We identify g%\f}J as the effective coupling between a wide class of excited 
string states on the pp-wave background. We compute the anomalous dimensions of 
BMN operators at first order in g\ and A' and interpret our result as the genus one mass 
renormalization of the corresponding string state. We postulate a relation between the 
three-string vertex function and the gauge theory three-point function and compare 
our proposal to string field theory. We utilize this proposal, together with quantum 
mechanical perturbation theory, to recompute the genus one energy shift of string 
states, and find precise agreement with our gauge theory computation. 
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1. Introduction 

Many years ago 't Hooft [1] demonstrated the existence of a nontrivial large N limit of 
SU (N) gauge theories 



In the 't Hooft limit (1.1), Yang-Mills interactions are controlled by the 't Hooft cou- 
pling A = g YM N. Away from the strict N — > oo limit, Yang-Mills perturbation theory 
may be organized as a double expansion. Feynman graphs are summed over their genus 
(controlled by the genus counting parameter 1/N 2 ) and over Feynman loops (controlled 
by the effective coupling A). These observations led 't Hooft to conjecture a duality be- 
tween large N gauge theories and weakly interacting string theories, 't Hooft proposed 
that the genus expansion on the two sides of this duality could be identified, leading to 
the identification of 1/N as the effective string coupling. The AdS /CFT conjecture and 
its generalizations have generated dramatic evidence for these proposals by supplying 
several concrete examples of such dualities. The study of these special examples has 
also led to the identification of A* as the effective string scale of the dual string theory, 
in units appropriate for comparison with the gauge theory. This implies, in particular, 
that as A — > oo, all string oscillator states have infinite mass and all unprotected single 
trace gauge theory operators have infinite dimension. 

Recently, Berenstein, Maldacena, and Nastase [2] have drawn attention to a differ- 
ent N — > oo limit of the N = 4, d = 4 Super Yang-Mills theory. The N = A theory 
has an 5*0(6) R symmetry group under which its six scalar fields X 1 . . . X 6 transform 



N 



°°> 9ym -> °> A = 9ym n fixed - 



(1.1) 
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in the vector representation. Consider an arbitrarily chosen U(l) subgroup of this R- 
symmetry group; for definiteness let this U(l) represent rotations in the X 5 and X 6 
plane. BMN study the sector of this theory with R charge J, and let J scale with N 
according to 



N 



• oo, with — ^= and Qym and A — J fixed. (1.2) 



Note that A — > oo and 1/N — > in the BMN limit. Consequently, according to 
the 't Hooftian lore reviewed above, SYM theory in the limit (1.2) is infinitely strongly 
coupled. Furthermore its string dual appears to be a free string theory with infinite 
effective string mass. None of these expectations is true; usual 't Hooftian reasoning 
fails as a consequence of the fact that observables in BMN limit are not held fixed, 
but scale to infinite charge as N — > oo. We will explain these remarks further below. 
However, it is useful to first review the string dual of Super Yang-Mills theory in the 
BMN limit. 

BMN were led to the large N scaling (1.2) by the consideration of a limit of the 
AdS/GFT duality. Super Yang-Mills in the seemingly singular regime (1.2) is actually 
dual to a well behaved closed string theory: IIB theory on the Ramond-Ramond pp- 
wave [4]: 

ds 2 = -4dx + dx~ - fi 2 z 2 dx +2 + dz 2 , -F+1234 = -F+5678 = — — ^, e* = g s . 

Aii A g s at z 

(1.3) 

According to this duality*, the R charge J of a Yang-Mills operator is proportional 
to the light-cone momentum p + of the corresponding string state, while A — J of the 
Yang-Mills operator is proportional to the light-cone energy p~ of the same state. The 
detailed dictionary between charges of the string theory and the gauge theory is given 
by 

/ip + a' = — =, = A - J, g\ M = 4ng s . (1.4) 

VA H 

Consequently, the AdS /CFT duality predicts that Super Yang-Mills theory in the limit 
(1.2) is dual to an interacting string theory with finite effective scale. This prediction 
is in conflict with the 't Hooftian expectations of the previous paragraph. 

We first address the puzzle of the effective string mass [2]. It is certainly true 
that all fixed unprotected single trace operators scale to infinite anomalous dimension 
(consequently the corresponding modes in the dual string theory scale to infinite mass) 

*This duality and its generalizations have been studied further by many authors, see [13]-[49]. 
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as A is taken to infinity. However, as we have emphasized above, observables are not 
held fixed, but scale with N in the BMN limit. While most such operators leave the 
spectrum in the N — > oo, A — > oo limit (1.2), BMN have identified a special set of 
operators whose anomalous dimension remains finite in this limit. These operators are 
dual to stringy oscillator states on the background (1.3). These operators are special; 
though they are not BPS, in the large N limit they are 'locally' chiral (see section two 
for more details), and so are nearly BPS. Scaling dimensions of these special operators 
do receive loop corrections; however the supersymmetric cancellations responsible for 
the non renormalization of exactly chiral operators also ensure that the anomalous 
dimensions of these almost BPS operators are much smaller than the power series 
in QymN that naive perturbative estimates suggest. Indeed BMN have argued that 
the anomalous dimensions of these special operators are not just finite, but actually 
computable perturbatively, even though the 't Hooft coupling A diverges in the limit 
(1.2). Supersymmetric cancellations produce a new coupling constant 

J 2 ~ (/ip+a') 2 ' { ] 

which appears to play the role of the loop counting parameter in the computation of 
two point functions of these operators. 

Like their scaling dimensions, three point functions of chiral operators are not 
renormalized [6, 7]. Consequently we expect analogous supersymmetric cancellations 
to permit the perturbative computation of three point couplings of BMN operators 
(hence interactions of the corresponding string modes) at small A'. In the rest of this 
paper (which is devoted to the study of PP-wave string interactions from perturbative 
Yang-Mills theory) we proceed on this assumption. The coherence and consistency of 
the picture that emerges provide some justification for this assumption. 

We now turn to the puzzle of the effective string coupling. String loops certainly 
contribute to scattering of modes of IIB theory at nonzero g s on the background (1.3) 
(see [8]), consequently generic correlation functions in Yang-Mills must also receive 
contributions from higher genus graphs even though N = oo, as in the limit (1.2). As 
we will demonstrate in section 3 of this paper, this puzzle has a simple resolution. It 
is certainly true that each graph at genus h is suppressed relative to a planar graph 
by the factor 1/N 2h . However we will demonstrate below that the number of diagrams 
at genus h is proportional to J Ah , so that the effective genus- counting parameter is 
actually 

^ 2 =(^) 2 = 16vr\ 2 (^W) 4 . (1.6) 
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g%, the effective genus counting parameter, must also control the mixing between single 
and multi trace operators; this is easy to see directly. The two point function between 
single trace and double trace operator is of order (fe/V^ (see section 3). A single trace 
operator of size J mixes with J different double trace operators; consequently this 
mixing contributes to two point functions at order J x ((fe/v^/) 2 = g\, in agreement 
with (1.6) for a genus one process. 

The identification of g\ with the Yang-Mills genus counting parameter fits natu- 
rally into duality between Yang-Mills and String theory, as gi has a rather natural 
interpretation in IIB theory on (1.3). In the pp-wave background, the worldsheet fields 
for the eight transverse directions are massive, so low energy excitations are confined to 
a distance 1/y/ np + from the origin. g\ = g 2 y/a'fip + is simply the effective two dimen- 
sional Newton's constant, obtained after a 'dimensional reduction' on the 8 transverse 
dimensions. 

In summary, despite first appearances, Yang-Mills theory in the limit (1.2) appears 
to develop a new perturbative parameter A'. In particular the theory is weakly coupled 
at small A' or large /i. Further, the genus expansion and mixing between single and 
multi trace operators — effects related to interactions in the string dual — are controlled 
by (?2, the effective two dimensional Newton's constant of the string theory. With 
this framework in place we proceed, in the rest of this introduction, to describe the 
precise relationship between string interactions and Yang-Mills correlators. As Yang- 
Mills correlators are perturbatively computable only at small A' or large //, some of the 
discussion that follows applies only to this limit. 

The first and most important qualitative issue concerns the identification of the 
effective string coupling in the background (1.3). Following our discussion of the genus 
expansion in gauge theory, it is tempting to identify the effective string coupling with 
g 2 . This guess is incorrect. In section 5 we will argue that the effective string coupling 
between states with the same A — J (at small A') in the pp-wave background is g2\/\' , 
where Ao is the scaling dimension at A' = 0. Note that the genus expansion of Yang- 
Mills theory (governed by the parameter g%) survives even in the free limit (A' = 0) 
when the effective string coupling is zero. This genus expansion appears to be rather 
unphysical; we believe it contains information about the map between string states 
and Yang-Mills operators, but does not appear to directly encode interesting stringy 
dynamics. Physical effects (like anomalous dimensions) from higher genus graphs are 
obtained only upon adding some Yang-Mills interaction vertices to these graphs; this 
addition leads to the re identification of the string coupling as s^v^A 7 - Note that string 
states in the pp-wave background blow up into giant gravitons [12] when J 2 /N ^> l/g s , 
i.e. precisely when g2y/X, the effective string coupling, is large. 
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Let us explain our proposal for 
string interactions in more detail. The 
spectrum of string states in the pp- 
wave background clumps into almost 
degenerate multiplets at large /i. The 
splittings between states of the same 
multiplet are of order fi\', while the en- 
ergy gap between distinct multiplets is 
of order \x. In section 5 we propose that 
the matrix element of the light-cone 
Hamiltonian between single and dou- 
ble string states within the same mul- 
tiplet is the three-point coefficient of 
the suitably normalized operator prod- 
uct coefficient of the corresponding op- 
erators (this quantity is O^g^)), mul- 
tiplied by the difference between their 
unperturbed light-cone energies. Since 
energy splittings within a multiplet are 
of order fi\ f , Hamiltonian matrix ele- 
ments between such states are of or- 



A-j= 3 ^ 



P jumps by fj^X 
perturbative in SYM 



IP jumps by /va 
nonperturbative in SYM 



A„J=1 



A„J=0 



Figure 1: Hierarchy of energy scales for large 
fj>p + a' . Only the transitions between states with 
the same Ao — J, i.e. with the same number of 
impurities, can be calculated perturbatively. The 
big jumps change the energy by an amount 1/A' 
times bigger and therefore they result from non- 
perturbative effects in Yang-Mills theory. 



der /i^A', corresponding to an invari- 
ant string coupling of order g2^/~X- For a class of BMN states we compute these matrix 
elements perturbatively in Yang-Mills theory. Note that transitions between states 
with different Ao — J involve large changes in energy; such transitions appear to be 
non-perturbative in the gauge theory (see figure 1). 

In the paragraphs above we have presented a specific proposal relating interac- 
tion amplitudes in string theory with correlation functions of the dual gauge theory. 
In the next two paragraphs we describe the evidence in support of our proposal. As 
we describe below, our proposal passes a rather nontrivial consistency check. Further 
we have also partially verified our proposal by direct comparison of three point func- 
tions (computed in Yang- Mills perturbation theory) with three string light-cone matrix 
elements (computed using light-cone string field theory). 

We first describe the consistency check on our proposal. In section 5 we compute 
the shift in dimension of a class of BMN operators to first order in A' and first order in g%, 
i.e. on the torus. We find that the anomalous dimensions receive non-zero corrections 
from the torus diagrams with a quartic interaction between "non-nearest neighbor" 
fields (see figure 10). The anomalous dimensions are proportional to g%\', the square 
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of the effective string coupling, and are interpreted as mass renormalizations of excited 
string states. We then proceed to recompute the mass renormalization of excited string 
states using second-order quantum mechanical perturbation theory. We obtain the 
light-cone matrix elements needed for this computation from correlators computed 
in perturbative gauge theory, utilizing our prescription described above. These two 
independent computations agree exactly, constituting a highly nontrivial "unitarity" 
check on the consistency of our proposals. 

In section five we also compare our proposal for string interactions with matrix 
elements of the light-cone Hamiltonian of string field theory. The light-cone Hamil- 
tonian is generated by a two-derivative prefactor acting on a delta functional overlap. 
In section 5 we demonstrate that the three-point function of three BMN operators, 
computed in free Yang-Mills theory, reproduces the delta functional overlap between 
three string states at large /i. We conjecture that, in the same limit, the prefactor of 
this delta functional overlap reproduces the second element of our formula for matrix 
elements (the difference between the unperturbed energies of the corresponding states). 
We sketch how string field theory predicts a modification of our prescription for the 
case of the operators involving D^Z and fermions. 

We conclude this introduction with a digression that may help to put our work in 
perspective. Yang-Mills/String theory dualities have hitherto been understood, even 
qualitatively, only in regimes of strong gauge theory coupling. For instance, it has long 
been suspected that confining gauge theories may be reformulated as string theories, 
with tubes of gauge theory flux constituting the dual string. However, as flux tubes 
emerge at distance scales larger than 1/Aqcd, their dynamics is nonperturbative in 
the gauge theory. More recently the Maldacena conjecture has established a duality 
between a conformal gauge theory (with a fixed line of couplings) and string theories 
on an AdS background. However these dualities are well understood only at large 
values of the gauge coupling. In this paper, utilizing the BMN duality, we have taken 
the first steps in explicitly reformulating an effectively weakly coupled gauge theory as 
an interacting sting theory (IIB theory on the pp-wave background at large /i). As 
perturbative gauge theories are under complete control, a detailed understanding of 
this extremely explicit duality holds the promise of significantly enhancing our under- 
standing of gauge-string dualities in general. 

The rest of this paper is organized as follows. Section 2 contains a review of 
subtle aspects of the BMN paper of importance to us. In section 3 we explain the 
counting that identifies g<i = J 2 /N as the genus counting parameter in free Yang-Mills. 
We also present the computation of planar three-point functions and torus two-point 
functions of BMN operators in free Yang-Mills theory. In section 4 we compute the 
torus contribution to the anomalous dimensions of BMN operators. In section 5 we 
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present our proposals relating Yang-Mills computations to amplitudes of the string 
Hamiltonian. We also present a nontrivial unitarity check of our proposals, and compare 
our proposals to string field theory. In section 6 we conclude with a discussion of our 
results and directions for future work. The reader who is uninterested in the details 
of perturbative computations of Yang-Mills correlators can skip from section 3.1 to 
section 5. In Appendix A we present a precise definition of a class of BMN operators. 
In Appendix B we prove that D-terms interactions do not contribute to the correlation 
function computations presented in this paper. In Appendix C we present a rigorous 
and self-contained derivation of two point functions of BMN operators. In Appendix 
D we present an alternative method for Yang-Mills computations. 

Note: As we were completing our manuscript, related papers appeared on the 
internet archive [9, 10, 11]. [9] overlaps with parts of sections three and four of our 
paper, while [10] overlaps with parts of section 3 and section 5.3 of this paper. Our 
results disagree with those of [9] and [10] in certain important respects. Unlike both 
of these papers we find non vanishing anomalous dimensions for BMN operators on 
the torus at first order in Yang-Mills coupling. We identify g2\^ rather than g 2 as 
the effective string coupling at large ji. As noted above, we have presented a rather 
non-trivial unitarity check of our proposals. We have also compared our proposal for 
the three-string vertex with the Green-Schwarz string field theory [8]. 



2. Preliminaries 

2.1 The BMN operators 

The simplest single-trace operator with R-charge J is 



J = J: -TrZ J , (2.1) 



where 



(2.2) 



This is a chiral primary operator, with scaling dimension exactly equal to J at all A'. 
According to the BMN proposal it corresponds to the light-cone ground state \0,p + ), 
where the map between parameters is given by (1.4). 

Other protected operators may be generated from J by acting on it with 5*0(6), 
conformal, or supersymmetry lowering operators. For example, by acting on J+1 with 
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a particular SO (6) lowering operator yields 



j 



1 



Tr (<PZ J ) , 



(2.3) 



Vn j+1 



where we have defined the complex combinations of the scalars: 



= 



X 1 + iX 2 
V2 



X 3 + iX 4 
V2 



(2.4) 



Oq is chiral with scaling dimension A = A = J + 1; it corresponds to the string state 
a o^|0;P + ); where aff = (clq + ia^)/\^2. To take another example, J+2 acted on by 
two distinct SO (6) lowering operators yields the protected operator 



which corresponds to the BPS string state Oq cZq |0,p + ). Proceeding in this manner, all 
protected operators (operators dual to supergravity modes) of the Yang-Mills theory 
may be obtained by acting on J , for some J, with the appropriate number of lowering 
operators of various sorts. 

As noted in the introduction, only protected operators remain in the spectrum as 
N is taken to infinity with gy M held fixed, in any sector of fixed charge J. However 
when J is taken to infinity together with N, it is possible to construct operators that 
are locally BPS. These operators consist of finite strings of fields (all of which are BPS) 
that are sewn together (in the trace) with varying phases into an operator of length 
J — > oo that is not precisely BPS. An example of such a near BPS operator is 



We will usually abbreviate this as J n \ however we must be careful to distinguish be- 
tween the two chiral operators Oq and Oq Q . In an inspired guess, BMN conjectured 
that the operator J n corresponds to the string state 0, p + ). As we have empha- 

sized above, for n ^ this operator is weakly non-chiral and its scaling dimension is 
corrected. However these corrections are finite, and may be expanded in a power series 
in A' (this result follows to low orders from direct computation, but independently, to 
all orders by comparison with the exactly known string spectrum). Operators corre- 
sponding to more than two string oscillators acting on the vacuum are discussed in 
appendix A. 




(2.5) 




(2.6) 
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J n was obtained from J+2 by replacing two Z's by the 'impurities' and ^, and 
sprinkling in position dependent phases. The impurities and ip were obtained by 
the action of 5*0(6) lowering operators on Z. In an analogous manner the impurity 
D^Z may be obtained by acting on Z with the generators of conformal invariance. 
Similarly, supersymmetry operators acting on Z produce gauginos. General BMN 
operators consist of these impurities sprinkled in a trace of Z's, together with phases. 
For the purposes of this paper it will be sufficient to consider only scalar impurities, 
but we will explain in section 5 how our ideas can be extended and checked with the 
other types of impurities. 

As we have stressed in the introduction, the dimensions of operators such as J n re- 
main finite (and perturbatively computable at small A') in the limit of infinite 't Hooft 
coupling only because these operators differ very slightly from protected chiral oper- 
ators. It is very important that the operator J n is defined to reduce precisely to the 
chiral operator Oq when n is set to zero. Even a small modification in the definition 
of this operator (such as a modification of the range of summation of the variable / 
to 1, . . . , J, as originally written in [2]) introduces a small — 0{l/\/~J) — projection onto 
operators that are far from chiral, resulting in perturbative contributions to scaling 
dimensions like gyyiN/J, which diverges in the BMN limit, and hence a breakdown of 
perturbation theory. ^ 

2.2 On the applicability of perturbation theory in the BMN limit 

Consider the perturbative computation of, say, the planar scaling dimension A of a 
BMN operator such as J n in (2.6) above. Suppressing all dependence on n, the results 
of a perturbative computation may be organized (under mild assumptions) as 

A = ±( 9 -^) m U9 2 YM N), (2.7) 

m=0 V J / 

where f m are unknown functions of the 't Hooft coupling. BMN computed the planar 
part of /i(0) using Yang-Mills perturbation theory, and deduced /i(oo) using the duality 
to string theory on the pp-wave background. Quite remarkably they found that fi (0) = 
/i(oo). This result suggests that f\(x) is a constant function at the planar level. 
Recently, the authors of [11] have demonstrated that /2(C)) = /^(oo), and have presented 
arguments which suggest that the planar components of f m (x) are constant functions 
for all m. Note that a term proportional to x m in f m (x) would result in the breakdown 
of perturbation theory, in the BMN limit, at order (g\ M N) m+n . Consequently, the 
conjecture that f m (x) are constant functions for all m is identical to the conjecture 

^The importance of the summation range 0, . . . , J has been also recognized by the authors of [9]. 
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that A' = 9ym^/ J 2 i s the true perturbation parameter, for the computation under 
consideration, in the BMN limit. 

In this paper we will proceed on the assumption that A' is indeed the perturbative 
parameter for the computations we perform, namely low order calculations of non- 
planar anomalous dimensions and three-point functions of BMN operators. We will 
see that (2.7) acquires extra non-planar contributions proportional to (J 4 /N 2 ) h from 
genus h diagrams. These contributions are finite in the BMN limit but they can be 
expanded in A' just like (2.7). All our results are consistent with this conjecture, and 
lend it further support; however, it would certainly be interesting to understand this 
issue better. 

3. Correlators in free Yang-Mills theory 

3.1 Correlators of chiral operators at arbitrary genus 

Consider a correlation function involv- 
ing operators of typical size (R-charge) J 
in free U(N) Yang-Mills theory.* Following 
BMN, we study this correlator in the large 
N limit; J is simultaneously scaled to infin- 
ity with J 2 /N held fixed. In this section we 
will demonstrate that the number of graphs 
that contribute to this correlation function 
at genus h scales with J like J Ah . Since any 
particular genus h graph is suppressed by a 
factor of 1/N 2h compared to a planar graph, 
we conclude that the net contribution of all 
Figure 2: Genus one diagram drawn on a S enus h graphs remains finite in the BMN 
square. limit, scaling like g% h where g 2 = J 2 /N. Con- 

sequently g\ is a genus counting parameter; 
it determines the relative importance of higher genus graphs in free Yang-Mills theory. § 
The free Yang-Mills genus expansion encodes a modification in the dictionary be- 
tween string states and Yang-Mills operators, but does not in itself appear to contain 
information about string interactions. We will return to the question of true string 
interactions in sections 4 and 5 below. 

Consider the two-point function (0 J (0)O J (x)) in free Yang-Mills theory, where the 

operators Q J are defined in (2.1). The planar contribution to this two-point function 

+Most formulae simplify for U(N) as compared to SU(N) and the relative difference is of order 
1/N. We therefore choose to work with the gauge group U(N). 

§It was previously observed in [52, 53] that operator mixing and higher genus contributions to 
correlation functions are important, even as N — * oo for operators whose size scales with N. [53] 
has also presented detailed formulae for free field correlators at all N in a basis different from that 
employed in this paper. 






Figure 3: The same genus one diagram in an alternative representation. The diagram can be 
interpreted so that (a) the blue string Tr Z 5+4+2+3 splits into Tr Z 5+ and TrZ 2+3 ; these two 
strings rotate by 5 or 2 units, respectively, to get the final state contracted with Tr Z 4+5+3+2 . 
The same free theory diagram however also counts (b) a similar split of the string into Tr Z +2 
and TrZ 3 + 5 . 



is 



(O J (0)O J (x)> 



planar 



(Ati 2 x 



2\J' 



(3-1) 



To find the genus 1 contribution to the correlator, we must find all the free diagrams 
that can be drawn on the torus but not on the sphere. To do this the J propagators 
must be divided into either 3 or 4 groups (see figure 2). The number of ways to do this 
is 



Mi 



j + i 

4 



J 4 
4! ' 



(3.2) 



This must be multiplied by J for overall cyclic permutations, but then divided by J 
again due to the normalization of the operator, and also by iV 2 due to the genus. The 
resulting quantity is finite in the BMN limit, and proportional to g\: 



<O J (0)O J (x)) toms 



9l 



4!(47r 2 a; 



2\J ' 



(3.3) 



This counting is easily extended to arbitrary genus. A genus h Feynman graph 
can be drawn on a 4/i-gon with sides identified pairwise. As we see from figure 4, the 
number of graphs that can be drawn on a 4/i-gon is the number of ways of dividing J 
lines into Ah groups, which is J 4h /(Ah)\. (The lines may also be divided into Ah — 1 
groups, but this gives a vanishing contribution in the BMN limit.) We must multiply 
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Figure 4: (a) Genus two surface represented as an octagon. The vertices of the octagon are 
all identified, while the edges are identified pairwise. The usual homology one-cycles with 
intersection numbers #(aj,6j) = 5ij are depicted. 

(b) Genus two diagram. One operator is located at the center, and the other at the vertices. 
This way of drawing the diagram immediately generalizes to any genus. 

this by the number of inequivalent ways of gluing the sides of a 4/i-gon into a genus h 
surface. This number has been computed [51]; the result is 

l'3--(4ft-0 (34) 
2h + l K J 

Consequently a total of 2~ 2h J 4h /(2h + 1)! graphs contribute to this correlator at genus 
h. Summing over genera we find 



(O J {0)O J (x)) 



(4tt 2 x 



2\J 



E 

h=0 



1 



(2/i+l)! V 2 



g 2 \^ 2sinh(<? 2 /2) 



g 2 (4:Tr 2 x 



2\J 



(3.5) 



This method can easily be generalized to show that the two-point function for an 
arbitrary chiral BMN operator such as Oq or Oq (defined in (2.3) and (2.6) respec- 
tively) has the same coefficient as in (3.5). Thus for example, 

2sinh((? 2 /2) 



<o yo)o yz)> 



(3.6) 



g 2 (4:ir 2 x 2 ) J+2 ' 

The easiest way to generalize to higher-point functions of chiral operators is prob- 
ably via a Gaussian matrix model. For example, it is not difficult to compute 



VZVZ 



Y[TiZ J ATiZ J e- Tl ^ 



,i=l 



j 



(3.7) 



i=l 



yielding simple explicit formulae that generalize (3.5).' 

^These formulae have been obtained in collaboration with M. van Raamsdonk. They have also 
been presented in detail in the recent paper [9]. 
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3.2 Planar three-point functions 

In this subsection we compute free planar three-point functions for the BMN operators 
defined in subsection 2.1. The results we obtain will be used in section 5 when we 
discuss the construction of string interactions. 
We will first compute 

(^(0)O^i)O j2 (:r 2 )>, (3.8) 

where J\ + J 2 — J ■ The planar, free field computation of this correlator is summarized 
in figure 12. The only complication is that we must sum over all of the possible 
positions for the and ip fields and carefully keep track of combinatorial factors as 
well as normalizations. The summation over the position of <fi and ip in O^ 1 may be 
converted into integrals in the large N limit, 

j\ C da f dbe ^Mm-ny) e -2nib(m-ny) = j2 ^ ™V (3 9) 

Jo Jo Tx\ny-my 
where y = J\/J. The final result for the correlator is obtained by multiplying this 
integral by J 2 (from cyclic rotations of Ji ) and dividing by y/ J x J 2 J (from the normal- 
ization of each operator) and by N (from 1/N counting). We find 

(O^^W) = r- . (3.10) 

A similar calculation yields 

(0^)0^)0^)) = — ^\nny) 

V nK ' ° 1 ^ ° 1 ;/ v / J7r 2 n 2 (47r 2 a;f)^+ 1 (47r 2 a ;|)^+i' 1 ; 

where Oq 1 and Oq 2 have and ip impurities respectively. 

These expressions for the three-point functions will play an important role in our 
comparison between perturbative string theory and perturbative Yang-Mills theory in 
section 5. 

3.3 Torus two-point functions of BMN operators 

In this subsection we present an explicit computation of the two-point functions 
for the BMN operators (2.6) at genus one in free Yang-Mills theory. The operator 
J n differs from the chiral operator Oq only in the presence of phases. Torus (and, 
indeed all genus) two point functions of Oq were computed rather easily in (3.6). The 
additional phases complicates matters somewhat, as we will see below. Nonetheless, it 
is not difficult to convince oneself that these additional complications affect only the 
details of the result, but not its scalings with J and N. Indeed, g\ is the correct genus 
counting parameter for all free Yang-Mills computations in the BMN limit. 
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I IL 



JL 



JL 



J 

J 5 



We consider first the correlator 
(O^(O)O^(x)). The free genus one diagrams 
are given by the torus diagrams presented in 
the last section (figure 3) with J + 2 lines, 
summed over all ways of replacing one line 
by a if) line and another by a line, with the 
rest becoming Z lines. There are four groups 
of lines, and if the if> and are in different 
groups then their relative positions will be 
different in the first and second operators, 
giving a non-trivial phase (unlike in the case 
of planar diagrams where they are always the 
same distance apart in the first and second 
operators). 

In fact it is convenient always to put the 
line at the beginning of both the first and 

second operator. With this convention every torus diagram with one line and J + 1 
other lines can be drawn as in figure 5, where each solid line represents a group of Z 
lines, and J\ + ■ • • + J 5 = J+l. Now we must put in the if> line. If it is in the first group 
(Ji possibilities) or the last group (J5 possibilities), then the phase associated with the 
diagram is 1, because it doesn't move relative to the line. On the other hand if it's 
in the second group ( J 2 possibilities) then it moves to the right by J 3 + J 4 steps, giving 
a phase exp(27rm(J 3 + J 4 )/J). Similarly for the third and fourth groups, giving in all 

Jt + J 2e ^in(J 3 +J 4 )/J + J 3e 2nin(J 4 -J 2 )/J + J 4 g-2™( J 2 + J 3 )/ J + Jg Wg mugt nQW gum ^ 

over all ways of dividing the J + l lines into five groups: 



Figure 5: A typical free genus one graph 
for the two-point function of an operator 
with two impurities. 



~ (^n(O)^n ( X ))froe torus 

92 

= -L J2 {Ji + J 2 e 2 ™ (j3+J4)/J + j 3 e 2 ™( J 4-J 2 )/J + Jie -2 m n(j 2 +j s )/j + 

J1+-+J5 
=j+i 



Af^oo 



Jo 



X (jl + J 2 e 2nin ^ + ^ + j 3e ^Mj4-j2) + j^ e -2nin(j 2 +j 3 ) + 



J_ 

24' 
1 



1 + 7 



. 60 6(27rn) 2 (2nn) 4 ' 



n = 0, 
n j£ 0. 



(3.12) 



In taking the limit — > 00 the fractions ji = Ji/ J go over to continuous variables. 
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If we now consider a correlator of two different operators J n and O^, then the 
phase associated with a diagram depends not just on which group the ip is inserted 
into, but on where in the group it is inserted. The formulae are therefore somewhat 
more complicated, but it's clear that again in the limit N — > oo the two-point function 
will reduce to g\ times a finite integral: 



(Aty 2 x 



2\J 



9l 



■ (0^(0)0 J m (x)) 

free torus 



N^oo 



1 djx ■ ■■dj b 5{j 1 H - 1) 



Au-v)h 



i(U-V) 

_|_ e iUji-iV(ji+j 3 +ji) ^ e i(U-V)j 2 _ 



1) 



J_ 

24' 

o, 

1 



6U 2 i [/4l 

4C7 2 13 ^ ' 

' (1 



' 60 



I (C/-V) 2 



4 4 _ D 

3 T y2 ~r JJ2 jjy 



_|_ e iU(ji+j2)-iV(ji+j4)^ e i(.U-V)j 3 
_|_ e iU(ji+j2+j3)-iVji ^ e i{U-V)j 4 _ 

+ e i(r/-y)0'i+j'2+j3+j4)^ e i(i/-v)j8 _ 
m = n = 0; 

m = 0,n^0orn = 0,m^0; 
m = n^0; (3.13) 
m = — n 7^ 0; 
all other cases 



(U-V) 2 



where U = 2-Km, V = 2irn. The result for the free two-point function including 
genus one corrections can thus be summarized as 



(O^(0)O J m (x)) 



free torus 



n 

(4ir 2 x 2 ) J 



(3.14) 



where the entries for A nm are given above. As (0^0^) is non-zero for n ^ m (unless 
either n or m is zero), we see that 3 n and J m mix with each other, and that the mixing 
matrix elements are 0(g 2 ). 

It is clear that the above procedure generalizes to the higher genus free diagrams 
described in section 2, in which the lines are divided into Ah groups. The genus h 
contribution to the two point function may be written as g^ times a finite integral 
over Ah + 1 parameters. See appendix C for a rigorous, general discussion. 



4. Anomalous dimensions from torus two-point functions 

The planar anomalous dimension of the operator J n is related, via the duality with 
string theory, to the light-cone energy (or dispersion relation) of the corresponding free 
string state. The planar anomalous dimension was computed to first order in g\ M in 
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[2]. Their result was of order O(X'), in precise agreement with the free spectrum of 
strings in the pp-wave background (1.3). On the other hand, the contribution to the 
anomalous dimensions from genus one gauge theory diagrams is related to the string 
one loop corrected dispersion relation for the corresponding state (see Section 5 for 
more details). In this section we compute the anomalous dimension of 3 n on the 
torus, to first order in g\ M . We find a result proportional to g%\', in accord with the 
identification of g\ as the gauge theory genus counting parameter, and A' as the effective 
gauge coupling. This result is a prediction for the one string loop 'mass renormalization' 
of the corresponding state. 

Below we present a diagrammatic com- 
putation of this anomalous dimension; in ap- 
pendices C and D two independent rigorous 
calculations confirm and generalize the re- 
sults of this section. 

We will find it convenient to think of 
the M = 4 Lagrangian in M = 1 lan- 
guage; Z, 0, ip are the lowest components of 
the three adjoint chiral superfields of this 
theory. Most of the interactions of the the- 
ory, including scalar-gluon (and ghost) inter- 
actions and scalar-scalar interactions of the 
form Tr | [Z, <p] | are 'flavor blind' (see Ap- 
pendix B). The contribution of these terms 
to this correlator is identical to their contri- 
bution to (Tr J J ); consequently they van- 
ish to order g Y M by the theorem proved in 
[7] (see Appendix B for more details). Con- 
sequently, only flavor sensitive terms in the Lagrangian, i.e. F-terms, contribute to our 
calculation. The F-term interactions between scalars are very simple 

V F = -4^ M Tr (p,0]| 2 + \[Z^)\ 2 + |[0,^]| 2 ) . (4.1) 

Further, at the order under consideration, the last term does not contribute, as it is the 
square of a term anti-symmetric under <-> ip and so has vanishing Wick contractions 
with 3 n and its conjugate, as these operators are symmetric in <f> and ip. In summary, 
to the order under consideration, the two impurities do not talk to each other, and may 
be dealt with individually. Further, each impurity effectively only interacts quartically 
with the Z fields through the interactions in (4.1). 




Figure 6: A typical graph with near- 
est neighbor interaction. Four (blue) lines 
coming from the vertex should be replaced 
by all possible terms from figure 7. The 
dashed line is a ip propagator. 
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Figure 7: The quartic Z-(j) vertices coming from the F-term (4.1). 

We now turn to the computation of all diagrams with a single F-term interaction. 
Consider contributions to the two point function 

(O J n (0)O J n (x)}. (4.2) 

from Feynman diagrams with a single Z, <p interaction vertex. All such graphs (see figure 
6 for one example) have identical spacetime dependence and their Feynman integral is 
proportional to 

1 f d*y _ ln(AV) 
16tt 4 J y\y - x) 4 8vr 2 a; 4 ' 1 ' 1 

We work in position space in (4.3); y represents the position of the interaction point, 
which must be integrated over all space. The integrand in (4.3) consists of two propa- 
gators from 0^(0) to the interaction point multiplied by two propagators from 0^(x) 
to the interaction point y, (see figures 6 and 2.1). 

In order to complete the computation of the torus two point function we must 

• (a) Enumerate all graphs that can be drawn with a single F-term interaction on 
the torus. 

• (b) Evaluate each of these graphs ignoring the propagators from the two operators 
to the interaction point (this corresponds to evaluating the corresponding free 
graph) and then multiply the result by (4.3). 

• (c) Sum over the contribution from all these graphs. 

In the rest of this section we carefully carry through this process to compute 
(O^(0)O^(x)) on the torus, to first order in the Yang-Mills coupling. It turns out 
that the graphs that contribute may be categorized into three separate groups; nearest 
neighbor, semi-nearest neighbor, and non-nearest neighbor graphs, respectively. 
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4.1 Nearest neighbor interactions 



Consider for example the diagram shown in figure 6, in which two adjacent lines in a 
free diagram such as figure 3 are brought together at an interaction vertex. We use 
the convention that diagrams at figure 6 actually represent the sum of four different 
Feynman diagrams. In diagrams such as figure 6, one of the lines connecting each 
of the operators to the interaction point is always <p propagator (two choices for each 
operator) and the other line always represents a Z propagator. The four Feynman 
diagrams correspond to the four possible choices. The dashed line on the right in figure 
6 represents a ip propagator. The four Feynman graphs that constitute the process 
depicted in figure 6 each contributes with the same weight; but graphs in which a 
line crosses the Z line contribute with a relative minus sign (this follows from the fact 
that the interaction is derived from g Y M Tr | [Z, 4>]\ 2 ), as shown in figure 7. The total 
contribution of these four diagrams is thus 



9ym | 
' N ' 



2nin/J 



)(1 



-2-7rm/ 



'n 2 



(2tt)V 



(4.4) 



times the phase associated to the correspond- 
ing free diagram. (4.4) is independent both 
of which two lines in the free diagram figure 
3 we are considering. It also does not de- 
pend on which particular free diagram is un- 
der consideration. Consequently, the sum of 
all such "nearest-neighbor" diagrams is sim- 
ply (4.4) multiplied by A nn , the genus one 
contribution to the free correlator (3.14) cal- 
culated in the previous section (with an ad- 
ditional factor of two from diagrams in which 
the interaction involves the ip rather than the 
field). 

Summing up all these diagrams, together 
with the free torus diagrams computed in 
this section, and adding these contributions 
to the free and one loop planar results com- 
puted in BMN we obtain the following cor- 
relator: 

(O J n (0)O J n (x)) = (47r ^ 2)J+2 (1 - A'n 2 ln(AV)) (l + g 2 2 A nn ) + 




Figure 8: A typical graph with semi- 
nearest neighbor interaction. 



(4.5) 
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Consequently, the diagrams studied in this subsection merely correct the coefficient 
of the logarithm in the two point function to account for the changed normalization 
of the operator as computed in the previous section. If there were no further 
contributions to the coefficient of the logarithm, this result would imply that torus 
diagrams do not contribute to the anomalous dimensions of the BMN operators." In fact 
other diagrams we describe in the next two sections do modify the scaling dimensions, 
as we describe below. 

4.2 Semi-nearest neighbor interactions 

There are two other classes of diagrams, 
illustrated in figures 8 and 9 (see also 
figure 10) that could potentially con- 
tribute to the correlator at order gfA', 
and thus to the anomalous dimension. As 
we will demonstrate below, the "semi- 
nearest-neighbor" diagrams of figure 8, 
in which the fields involved in the in- 
teraction are adjacent in one but not 
the other operator, contribute to the two 
point function (O^(O)O^(x)) only when 
m n. However, the "non-nearest neigh- 
bor" diagrams of figure 9 contribute to 
the logarithmic divergence of this correla- 
tor whether or not n = m. Consequently, 
these diagrams result in a genuine shift in 
the anomalous dimension of J n . 

It is not difficult to argue that no other classes of diagrams contribute to this 
process. To verify this claim, consider all diagrams with a single quartic interaction, 
that can be drawn on a torus. Each such diagram involves two propagator loops 
involving the interaction point. All diagrams fall into four classes; diagrams in which 
each of these loops is contractible, in which one loop is contractible and the other wraps 
a cycle of the torus, in which both loops wrap the same cycle of the torus, and finally 
those in which the two loops wrap different cycles of the torus. Further dressing these 
diagrams with all sets of propagators that leave it genus one, we find that the first class 
constitutes nearest neighbor graphs of the form figure 6, the second set constitutes 
semi- nearest neighbor graphs of the form figure 8, the third set constitutes nonnearest 
neighbor graphs of the form 9 and the last set cannot be implemented with F-term 
interactions. In Appendix C and D we verify this result using different techniques. 
" The unphysical nature of these contributions to the two-point function was also recognized in [9] . 




Figure 9: A typical graph with non- nearest 
neighbor interaction. Those graphs contribute 
to the actual anomalous dimension. 
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In this subsection we discuss the semi-nearest-neighbor diagrams. There are exactly 
eight diagrams of this type corresponding to the number of ways one may choose the last 
member of a given group in figure 8 to interact with the first member of the next group. 
The number of semi-nearest neighbor diagrams is smaller than the number of nearest 
neighbor diagrams by a factor of 0(1/ J) as either or ip must be located at the edge 
of one of the four 'groups' of lines in figure 8. Consequently such diagrams are naively 
negligible in the J — * oo limit. However, each individual semi- nearest neighbor diagram 
is enhanced by O(J) relative to a nearest neighbor diagram. In order to understand 
this, consider for example the case illustrated in figure 8. As in the previous case, this 
figure really represents four diagrams, which contribute a total 



(the last factor is due to the ip field, which in this particular example happens to sit in 
the fourth block, but whose position should be summed over). The fact that there is 
only one power of J in the denominator, rather than two as in (4.4), compensates the 
fact that these diagrams are rarer by a factor of 1/J than the nearest-neighbor ones. 
Consequently, such diagrams could make non-vanishing contribution in the BMN limit 
J — > oo, and they do contribute to (O^(O)O^(x)) for m ^ n. However it turns out 
that the full contribution from semi-nearest neighbor graphs to the correlator above 
vanishes for the case m = n considered in this section. One can see this either by 
considering the other semi-nearest neighbor diagrams at fixed J±, J 2 , J3, J4 (there are 
32 such diagrams in total), and seeing the cancellations explicitly, or by considering a 
fixed diagram such as figure 8 summed over J\, J2, J3, J a — since (4.6) is antisymmetric 
under exchange of J\ and J 2 it must vanish in the sum. 

4.3 Non-nearest neighbor interactions 

Finally, we turn to the non-nearest-neighbor graphs described in the figure 9 (redrawn 
differently in figure 10). The external legs of our operator are divided into three groups 
containing J x and J 2 or J 3 Z's, respectively. Because we have divided the Z propagators 
into three rather than four lines, these diagrams are rarer still by a factor 1/J than the 
semi-nearest-neighbor diagrams, or by a factor of 0(1/ J 2 ) compared to nearest neigh- 
bor diagrams. However, this is compensated by the fact that each non-nearest neighbor 
diagram is enhanced by a factor O(J) compared to semi-nearest neighbor diagrams, 
or 0(J 2 ) compared to non- nearest neighbor diagrams. In order to see this note that 
in the diagram of figure 9 both ends of the <fi propagator jump by a macroscopic (i.e. 
order J) distance along the string of Z's. Let the impurity be located on the left 
and/or right end of the group J\ whose first and last propagators are "pinched". The 




,2 



(4.6) 
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four diagrams represented by figure 9 consequently contribute equally, but weighted by 
phase 1, (for the two diagrams in which does not jump) or — exp(+2ninJi/ J) for 
the diagrams in which <fi jumps either to the left/ right; consequently the sum of these 
four diagrams is proportional to 



[1 - e 2ninJl/J )(l - e 



2mnJ\j J\ 



We now turn to the contribution to the 
diagram from the phase associated with the 
second impurity ip. If ip is in one of J\ 
places inside the first vertical block its rel- 
ative position on the two operators is the 
same, and so these J\ diagrams contribute 
with no phase. On the other hand, if if) is 
in the block with J 2 propagators; its relative 
position on the two operators is different by 
J3; the corresponding J 2 diagrams contribute 
with phase exp (2-7^71/3/ J). Finally, if ip can 
be located in the third block (with J 3 prop- 
agators) its relative position on the two op- 
erators slides to the left by J 2 units. Conse- 
quently, these J3 diagrams are each propor- 
tional to exp(— 2mnJ 2 / J). Replacing the sum over J 1 (with J\ + J 2 + J3 = J) by an 
integral over j% = -hi J (with j\ +32 + 33 — 1)> we arrive at the integral 




Figure 10: Figure 9 represented as a pe- 
riodic square. 



djxdhdhKh +32 



1 5 

3 + 2n 2 n 2 



(for n 7^ 0). The two-point function is thus, at first order in A' and g%, given by 



(4.7) 



(^x 2 ) J+2 (d J n (0)O J n (x)) 



1 + 9lA nn ) (1 - n 2 A' ln(AV)) + M Q + J-] l n(A V), (4.8) 



2n 2 n 2 J 



and the anomalous dimension is given by 



A = J + 2 + X'n 2 



gi* n 



4tt 2 V3 2n 2 n 2 



(4.9) 
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The methods described in this section, or the ones described in appendices C and 
D, may also be used to calculate the two-point function between different operators. 
We summarize the result here; the reader will find the details in the appendices: 

(5 nm + g\A nm ) (f - (n 2 - nm + m 2 )\' ln(AV)) + ^B nm ln(AV). (4.10) 

Here the first, factorized term contains the contributions of the nearest-neighbor (pro- 
portional to nm) and semi- nearest-neighbor (proportional to (n — m) 2 ) diagrams. The 
second term contains the contribution of the non-nearest-neighbor diagrams: 

0, n = or m = 0; 

5 + $» n = m^0; ^ 

n = —m 7^ 0; 



Burn 



15 

2C7 2 "' 



^ + (u_ V )2 , all other cases, 

where U = 2nm and V = 2nn. 

The classification of diagrams into nearest-, non- nearest-, and semi-nearest- 
neighbor continues to be valid at higher genus (at first order in A'). Interestingly, 
the factorization of the first two contributions, as in (4.10), is true at all genera. 



5. String interactions from Yang-Mills correlators 

In this section we finally turn to the relationship between correlation functions in Yang- 
Mills and dual string interactions. We make two specific proposals at large //: 

• Three-point functions of suitably normalized BMN operators, multiplied by the 
difference in p~ between the ingoing and outgoing operators, may be identified 
with the matrix elements of the light-cone Hamiltonian between the corresponding 
one string and two string states. 

• The one string loop mass renormalization of a class of excited string states is 
reproduced by the 0(g 2 ) anomalous dimensions of the corresponding operators. 

We believe that these proposals form part of a larger dictionary relating the Yang-Mills 
theory and the string theory; however we leave the determination of the rest of this 
dictionary to future work. 

This section is organized as follows. In subsection 5.1 we motivate and explain 
our proposals in detail, and elaborate on some of their consequences. In the rest of 
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this section we provide evidence for the validity of our proposals. In subsection 5.2 we 
demonstrate that our proposals pass a nontrivial self-consistency check. In subsection 
5.3 we compare our proposals (together with the computations of Yang-Mills correlators 
in section 3 and 4) with the predictions of string field theory, and find substantial 
agreement. 

5.1 Three-string light-cone interactions from Yang-Mills three-point func- 
tions 

Let Oi, Oj, and Ok represent three single-trace BMN operators, of U(l) charges Jj, Jj, 
and Jfc, and normalized so that 

(diW S (x)) = (5.1) 

Let represent the free single string states that correspond to these operators at zero 
bulk string coupling, normalized such that 

W) = (5-2) 

Let 

oji,jj+j k Cijk 



(O i (a; i )O i (a; i )C' fc (a; fc ))f re e p i a 



(27r^ j ) Al+Aj ~ Afc (27r^fc) Al+Afc ~ Aj (27ra; ifc ) A J + A ^- A ^ 

(5.3) 



The coefficients have been evaluated in (3.10) and (3.11). At small A', we pro- 
pose the following formula for the matrix element of the string field theory light-cone 
Hamiltonian 

(i'\p-\j'k') = //(Aj - A j - A k )C ijk . (5.4) 

(5.4) is expected only to apply to leading order in A'; we leave its generalization to 
finite A' to future work. 

Equation (5.4) is one of the central proposals of our paper. In sections 5.2 and 5.3 
below we will provide rather strong evidence for its validity. Before proceeding to do so, 
however, we provide initial motivation for the proposal (5.4). Inner products of Yang- 
Mills states on S 3 (and so, presumably, states of the dual string theory) are related 
to correlation functions of the Euclidean Yang-Mills theory by the state operator map. 
Thus it is plausible that matrix elements of the string theory light-cone Hamiltonian 
are given by Yang-Mills correlators, dressed by a factor of linear homogeneity in p~ . 

We now motivate the specific form of the dressing in (5.4). Yang-Mills correlators, 
correctly normalized (see below), are of order g\. On the other hand, from section 4, 
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torus mass renormalizations occur at order /zgf A', and so go to zero when A' is taken 
to zero at fixed g 2 - Consequently, the dressing factor must go to zero as A' is taken to 
zero; this suggests the specific form of the formula (5.4). 

In the rest of this subsection we elaborate on the consequences of (5.4). 

5.1.1 Scaling with N and J 

Note that C^k scales like J 3 ^ 2 /N for the BMN operators under consideration.** Further, 
in the large /i limit the energy splittings (p± + p 2 — p% ) are of order //A'. Consequently, 
the right-hand side of (5.4) scales like /i^A' / \fj. As J is taken to infinity in the BMN 
limit, these matrix elements scale to zero, which is puzzling at first sight. Note, however, 
that the number of intermediate states (or final states) in any process scales like J; 
consequently (see subsection 5.2) the scaling of matrix elements is precisely correct to 
yield finite contributions to physical processes. Stated differently, the scaling of matrix 
elements like //g^A' / \fl is merely a consequence of dealing with string states that are 
unit normalized. Switching to the more conventional delta function normalization for 
states 

m=pts(pt-pt) = j l s. h ,j J (s.e) 

requires a rescaling of states 

\i) = y/Ttf)- (5.7) 
Light-cone Hamiltonian matrix elements may then be written as 



(i\P~\jk) 



(P, ~Pj -Pk)\\ J - 2 J LC nk 



pt s (pf -pt -pi) ( 5 -* 



The term in the square bracket on the RHS of (5.8) is finite in the BMN limit and is 
of order /zg^A' '. 



**This is easiest to verify in a simple example. The normalized chiral operators O j = Tr Z 1 ' j\J N J J 
have planar three-point functions 

<O J (0)O Jl (*i)O Ja te)W - (^P^gfcrj (5-5) 
where C Ju j 2tJ = y/JiJ 2 J/N. 
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5.1.2 Effective Coupling 

It is instructive to perform the following exercise. Consider an effective two dimensional 
field theory with scalar fields fa, interacting through a p + dependent cubic interaction 

oc J dx + dpfdp+dp+g(pf,p+,p+)5(pt+p+ - p+ )fa(x + ,pt)(f) j (x + ,p^)(j) k (x + ,p'l) 

(5.9) 

Canonically quantizing this theory in the light-cone, it is not difficult to verify (for ex- 
ample, by adapting equation (23) of [50] to our normalization) that the matrix elements 
for the light-cone Hamiltonian of this system are 

(i\P~\jk) oc J dpjdp+dplg(pi,p+,pl)5(pj +p+ -pi) (5.10) 

Consequently we conclude that (5.8) would be reproduced from a two dimensional cubic 
effective field theory with coupling (of dimension squared mass) given by g(pf, Pj~,p%) — 
(Pi ~Pj -Pk)V J j J k/ J i C ijkPt, i-e. 



9(pt,pt,Pt) - V 1 + A 2 - A3)^Mc 123 ~ (5-11) 

leading to the identification of (?2V / A' as the effective string coupling for these processes. 
5.1.3 Vanishing of on-shell amplitudes 

Recall that, in field theory, the decay of a particle is the result of the mixing between 
single particle states and multi particle states of the same energy. This mixing invali- 
dates the use of non-degenerate perturbation theory in following the 'evolution' of the 
unperturbed single particle state upon turning on an interaction. It fuzzes out the 
very notion of a particle; in particular the mixing broadens out delta function peaks in 
spectral functions, endowing the 'particle' with a finite lifetime. 

It is striking that (5.4) prescribes the vanishing of matrix elements of the light-cone 
Hamiltonian between states of equal unperturbed energy. This prescription implies the 
stability of excited string states in the large \x limit even upon turning on interactions. 
As the notion of a single particle continues to be well defined in the interacting the- 
ory, it is thus natural to identify the BMN operator (2.6) with the stable one-particle 
state, at large \i, even upon turning on interactions.^ This feature (the vanishing 
of matrix elements between states of equal unperturbed energy) also permits the use 



ttit may be possible to derive this identification, together with our proposal (5.4), from a careful 
analysis of the state-operator map. We hope to return to these issues in the future. 
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of non-degenerate quantum mechanical perturbation theory in an analysis of mass 
renormalization of excited string states at large //. We will utilize this observation in 
subsection 5.2 below. 

5.2 Unitarity check 

As we argued in the introduction, eight transverse coordinates in the pp-wave back- 
ground are effectively compactified. The light-like direction is also compact at finite 
N (its conjugate momentum, J, is quantized) and string theory on the pp-wave back- 
ground reduces to quantum mechanics. In this subsection, we apply standard quantum 
mechanical second order perturbation theory to perform a self-consistency check of the 
amplitudes calculated from the gauge theory. The Hamiltonian here is A = J + P~//i. 

Consider the string state corresponding to the BMN operator O'l defined in equa- 
tion (2.6). We will use the well known formula for non-degenerate second order per- 
turbation theory 

^ = E <« 2 > 

to compute its second order energy shift. 

In (5.12), the sum over states m includes two types of intermediate states: 

• (A) the two-string states with strings corresponding to and O j2 (these must 
be summed over the worldsheet momentum m). Using (5.4) and (3.10), the 
squared matrix element that connects 3 n to this two-particle state is 

, T/ i2 _ d^i 1 ~ y)( n V + m f sin 4 (7rm/) 

\V mn \ — 4 T ( (O.lOj 

where we have defined y = J\j J . The difference in energies between our state 
and the two-particle state is 

E n -E m = X > V - m2 \ (5.14) 

y 

• (B) the two-string states described by the two chiral primaries Oq 1 and Oq 2 , 
where the impurities are and ip in the two operators respectively. The squared 
matrix element of the light-cone Hamiltonian connecting J n to this two-particle 
state is easily computed from (5.4) and (3.11): 



IT J 



The difference in unperturbed energies between the initial and intermediate states 



is E n = X'n 2 . 
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In both cases, we must sum over J\ = J — J 2 , he. integrate over y — J\j J from to 1. 
The total torus correction to the dimension of 3 n is therefore 

Ai 2 ) = j/ 1 dyfsW + X;^ ) V (5-16) 



where 



A (n) _ \Vmn\ 2 = gl^'yil - y)(ny + m) sin 4 (7rm/) 
m ~ E n -E m " './(//// - ///)•' 



and 



(5.17) 



D (n) _ iKol 2 4g|VsinSmy 

5 = = ttv^ (5 - 18) 

The sum over m in (5.16) may be performed using the identity 

ETiy -\- Tfl 22/\/ / \\ / \ 
= —71 csc limy) (1 + 2imy coti-nny)) . (5.19) 
(ny — m) 6 

Adding to the result and performing the integral over y, we find 

^ = "Sr g + ^) • <**>) 

which is precisely (4.9), the genus one contribution to the anomalous dimension. 

In conclusion, we have computed the one loop mass renormalization from gauge 
theory in two different ways; from the genus one contribution to the anomalous di- 
mension of the corresponding operator, and independently using our proposal (5.4) 
and standard second order perturbation theory. These two computations agree ex- 
actly. In the next subsection we will proceed to compare our prescription (5.4) with 
the predictions of string field theory. 



5.3 Comparison with string field theory 
5.3.1 The delta-functional overlap 

In this subsection we show that the free planar three-point function in Yang-Mills 
theory is identical to the delta-function overlap between string states in the large \i 
limit. 



- 28 - 




Figure 11: A Feynman diagram for the planar three-string correlator and the corresponding 
light-cone gauge history of joining strings. 

In the light-cone gauge, the bosonic part of the worldsheet action of a string prop- 
agating in the pp-wave background is 

(d t X l d t X l - deX'drX* - t i 2 X i X i ) . (5.21) 

In the limit A' — > 0, at finite n the sigma derivative piece in (5.21) (d cr X t d cr X l ) is 
negligible compared to the mass term, and may simply be ignored, to zeroth approxi- 
mation. It is convenient to discretize the worldsheet (5.21) into J bits. On neglecting 
d a X l d a X t , the bits decouple from each other, and the string disintegrates into J inde- 
pendent harmonic oscillators. 

Now consider an interaction process involving three of these strings. String split- 
ting/joining interactions in light-cone gauge are local, and the most important piece of 
the three-string interaction Hamiltonian is the overlap delta functional 

J VX\a)VX\a)VX\a) A(X l (a) + X 2 (a) - X 3 ( ( x))^(X 1 )^(X 2 )^(X 3 ). (5.22) 

On discretizing the three strings, every bit on the third string (the biggest of the 
three) is put in correspondence with a bit on one of the two smaller strings. Specifically, 
two bits are in correspondence if they share the same value of a. The wave functions 
ip{X) factor into wave functions for each bit. Any bit is either in its harmonic oscillator 
ground state (we could denote that, in a diagram, by writing the letter Z in the 
appropriate slot) or in the first excited state of i th harmonic oscillator (denoted by 0* 
in the appropriate slot). The interaction (5.22) is simply 

V = H{x'i\Xi), (5.23) 
i 

where the index I runs over the bits of the larger string, is the harmonic oscillator 
wave function of the I th bit on the larger string, and \x[) is the harmonic oscillator wave 
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function on the corresponding bit on the corresponding smaller string (either string one 
or string two, depending on the value of a). Diagrammatically, V is unity when Z's sit 
on top of Z's and </>*s on top of 4> l s. It is zero otherwise (see figure 12). 



This is precisely the rule to compute the 
free planar contribution to the three-point 
function term by term in the series for the 
operators (2.5). The sum over phases in 
(2.5) is just the discrete Fourier series which 
may be taken simultaneously in Yang-Mills 
theory and on the discretized string world 
sheet. It is thus guaranteed that the results 
(3.10) and (3.11) from free planar gauge the- 
ory diagrams precisely reproduce the delta- 




functional overlap of strings in the /i — ► oo 

limit. Figure 12: Three-point function in free 

theory as the delta functional overlap of 
5.3.2 The prefactor the initial and final state. 

Our prescription (5.4) for the matrix elements of the light-cone Hamiltonian involved 
two elements. 

• (a) The free Yang-Mills correlator. 

• (b) The dressing factor (p± — pT, — p%). 

Similarly, the one/two string light-cone Hamiltonian involves two elements, a delta 
function overlap, and a prefactor acting on that delta function overlap. In the previous 
subsection we have demonstrated that the delta function overlap is precisely equal to 
the free Yang-Mills correlator. In this subsection we conjecture that the action of the 
prefactor on this delta function precisely produces the dressing factor (p^ — p^ — pj) 
multiplying the delta function overlap in (5.4). 

The prefactor of string field theory involves two derivatives, and might naively have 
been estimated to be of order p + jj, at large \i. In fact, as \i — > oo, the prefactor vanishes 
when acting on the delta functional overlap!** This striking result is consistent with 
our proposal that the prefactor, acting on the delta functional produces a factor of 

PI-P2 -Pa =0(l/fM) (5.24) 

An honest verification of (5.24) appears to be algebraically involved, though it is 
straightforward in principle. 

■f+This result has been derived by M. Spradlin and A. Volovich, and will be presented elsewhere. 
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A more detailed understanding of the structure of the string field theory prefac- 
tor in the large (j, limit would permit several predictions for gauge theory correlators. 
For example, a preliminary analysis appears to indicate that the normalization of the 
prefactor (5.24) is proportional to the number of scalar impurities minus the number of 
D^Z impurities (with no contributions from fermionic impurities), and would seem to 
imply that the operators with one scalar and one D^Z excitation, as well as operators 
with fermionic excitations only, should have vanishing amplitudes. 

6. Conclusions and outlook 

In this paper we have begun an investigation into the relationship between Yang-Mills 
correlators in the BMN limit and string interactions in the pp-wave background. Our 
analysis is valid at large a'p + fi ~ (A')~ 1/2 where the gauge theory appears effectively 
weakly coupled. We have employed perturbative Yang-Mills theory to make verifiable 
predictions for interaction amplitudes and mass renormalizations of weakly coupled 
strings on the pp-wave background. 

The principal observations and proposals of our paper are 

• For the appropriate class of questions, it appears that Yang-Mills perturbation 
theory in the BMN limit may be organized as a double expansion in an effective 
loop counting parameter A' = gy M N/ J 2 [2] and an effective genus counting pa- 
rameter g\ = J 4 /N 2 . In particular, graphs of all genus continue to contribute, 
even in the strict M — > oo BMN limit. The effective genus counting parameter 
g\ may independently be identified as the two dimensional Newton's constant for 
the dual string theory. 

• Mixing effects between single and multi trace Yang-Mills operators are also con- 
trolled by the parameter g^. This implies a modification of the dictionary between 
Yang-Mills operators and perturbative string states, proposed in [2] at the same 
order. 

• We have proposed a relationship between Yang-Mills three point functions and 
three string interactions at large \i. Our proposal, (5.4), implies that Light-cone 
Hamiltonian matrix elements between a wide class of single and double string 
states are of order /x^A', corresponding to invariant effective string coupling of 
order g2\f\' ■ At large values of this coupling string perturbation theory breaks 
down and string states blow up into giant gravitons. The detailed form of (5.4) 
also implies that, at large /i, excited string states are stable, even at nonzero 
values of the effective string coupling. 



-31 - 



• We have computed the one loop correction to the dispersion relation of these 
string states in two different ways: first by a direct computation of the anomalous 
dimension associated with this operator at order g\ (i.e. on the torus), and second 
from quantum mechanical perturbation theory, using matrix elements obtained 
from three-point functions, as in the previous proposal. These computations agree 
exactly; this constitutes a non-trivial check on our proposals. They also confirm 
our identification of g%\' as the effective theoretic genus counting parameter. 

This paper suggests several directions for future investigation. To begin with, it 
would be useful to check the proposals of this paper, and to better understand some 
its assumptions. It is very important to check our proposal for the dictionary between 
Yang-Mills correlators and string interactions; to this end an explicit expression for the 
first term in an expansion in powers of of the string field theory interaction vertex 
formally derived in [8] is required. It is certainly important to thoroughly understand 
when and why perturbative Yang-Mills can be employed in the study of this strongly 
interacting gauge theory (see [11]). Finally, in the unitarity check of section 5.2 we did 
not account for intermediate states with different numbers of impurities from those in 
the initial state. Contributions from such intermediate states are suppressed by large 
energy denominators. However, they could also be enhanced by parametrically large 
couplings. It would be useful to understand precisely when and why it is justified to 
omit such contributions. 

Several generalizations of our work immediately suggest themselves. It should be 
straightforward to generalize our calculations and proposals to BMN operators involv- 
ing D^Z and fermionic impurities. More generally, it would be very interesting to 
extend the dictionary between Yang-Mills and string theory proposed in this paper, 
beyond the large \x limit, and to wider classes of observables. It may also be possible 
to derive (5.4) and its generalizations from the more usual AdS/CYT prescriptions. 

In conclusion, the BMN limit appears to allow us to re-interpret a sector of an effec- 
tively perturbative Yang-Mills theory as an interacting string theory! This remarkable 
idea certainly merits further study. 
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A. Specification of operators 

The operators of the M = 4 SYM theory which map into low-energy string states in the 
pp-wave background involve a large number J of scalar fields Z together with several 
impurity fields <fi and ip. It simplifies things to take these fields to be holomorphic 
combinations of the six real elementary scalars X % of the theory, e.g. Z = (X 5 + 
iX 6 )/V2, = (X 1 + iX 2 )/V2, ijj = (X 3 + iX 4 )/V2. We take Z, 0, ip to be the lowest 
states of the three chiral multiplets $ 1 , $ 2 , $ 3 which appear in the M = 1 description of 
the M = 4 theory. BMN have outlined the rules for the construction of these operators, 
but there are some subtleties, and we thus briefly describe a clear specification here. 
The principles of the construction are the following: 

• (1) For the case of / impurity fields there are / initially independent phases 
qj = exp(2ninj / J). One writes a formal definition of the operators which reduces 
to a BPS operator when all qi = 1, specifically a "level /" SU(4) descendent of 
the chiral primary operator Tr(Z J+I ). Many of the operators so defined vanish; 
specifically they vanish unless the product qiq2---qi = 1 which is the level- 
matching condition on string states. 

• (2) On the "level-matching shell" the non-vanishing operators still satisfy the 
BPS property when all qi = 1. They reduce to the familiar BMN operator 
(with a minor but necessary change) for k = 2 and extend the construction to 
general k. Planar diagrams for the 2-point correlation functions of these operators 
vanish (both for free fields and order g\ M interactions) unless the momenta rij 
are conserved. This planar orthogonality property is approximate, holding to 
accuracy 1/J in the BMN limit N — > oo, J ~ \^N. Since the momenta do not 
have any clear meaning in the field theory, one should not expect strict momentum 
conservation. Indeed non-conservation becomes a leading effect for diagrams of 
genus > 1, again both for order g\ M and order g\ M . 
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We start the discussion with the 1 = 2 BMN operator 

j 

Oj = Y t Q lr ^(^Z J - 1 ), (A.l) 

1=0 

modified so that the sum begins at / = rather than / = 1 as in [2]. With this definition 
the BPS property is exact when q = 1. This minor change is significant when one 
computes the planar correlator including interactions. For the operator defined above 
one finds by techniques described elsewhere in this paper that 

j 

(Oj(x)Oj(O)) = fj(x)[l + cg 2 YM N(l -q)(l- f) In \x\] J>f)<, (A.2) 

1=0 

where q = exp(2nin/ J) and r = exp(2nin'/J) are the phases of the operators Oj(x) 
and 0/(0), respectively, and fj(x) = (N/4ir 2 x 2 ) J+2 involves the product of free scalar 
propagators. The factor (1 — q)(l — f) is the discrete second derivative of the phases 
coming from the stepping effect of the interactions noted in [2]. This factor contributes 
a 1/J 2 suppression in the BMN limit. With the sum in the operator beginning at I — 1 
one would obtain a similar expression with the changes: a) the final sum in (A.2) starts 
at I — 1, and b) there is an additional term proportional to gy M N{l + qf) In The 
last term arises because the construction of discrete second derivative is incomplete at 
one site. It is not suppressed in the BMN limit. If present the physical interpretation 
would be spoiled, which is why is we chose the definition (A.l). 

The final sum in (A.l) has the value J+l if qf — 1 and the value 1 otherwise. In the 
BMN limit (to accuracy 1/J) one thus obtains the momentum conserving result JS n y 
as expected for a string in the state a*(n)b*(— n)\0) in the state dual to the operator 
Oj in the limit (/ip + — > oo) when string interactions vanish. 

The operator for two identical impurities is obtained simply by replacing ip by 
<f> in (A.l). There are additional cross terms in the diagrams for the 2-point function. 
The planar result is just (A.2) with the final sum replaced by 

j 

W + (or) 1 ] = J (5 ny + 5„,_n0 (A.3) 

1=0 

in the BMN limit. This is the correct orthogonality condition for (non-interacting) 
strings in the level-matched state a*(n)a*(—n)\0). 

We now turn to the general program of defining operators off the "level-matching 
shell" . The general method is embodied in figure 13 which shows a circular array of J 
points, the J Z's in the trace, with and ip in interstitial positions at distances k and 
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I from an arbitrary origin. The associated phase is q\q\. It is clear from the figure that 
if we fix the relative distance, say / between if) and 0, and sum over rigid displacements 
of the positions of from k — to k — J— 1, we accumulate the phase polynomial 
q l 2 [l+qiq2+- ■ ■ (<?i<?2) J_1 ] which vanishes unless the level-matching condition rii + n 2 = 
is obeyed. 

The analytic expression for the op- 
erator depicted in figure 13 involves a 
sum over the two relative orders of 0, ip 
within the trace. It is 

= 0! + 2 ; 

o 1 = q 1 tq l 2 r ^{z k 4>z l - k ijjZ J - 1 ); 

0<fc<«<J-l 

2 = qi4MZ l 'tpZ k '- l '<j)Z{A% 

0<l'<k'<J-l 

It is straightforward although a bit 
awkward to carry out analytically the 
sum over rigid displacements of a con- 
figuration of fixed relative phase I pic- 
tured above. Cyclicity of the sum is 
vital, of course. Starting from k = 
in Oi of (A. 4) and moving to position 
k = J— 1 — 1, one accumulates the phase 

polynomial q 2 [l + qiq2~\ ("Ji^)* 7- ' -1 ]- 

The next step takes us into 2 with 
/' = 0, k' = J — I, and we sum over 
I' in I — 1 rigid steps until <fi is in next-to-last position in the trace. The phase poly- 
nomial from this traversal of 2 is q(~ l [l + q\q% + • • • (gi^)' -1 ]- The sum of these two 
polynomials is equal to the full polynomial in the previous paragraph, thus giving the 
level- matching condition exactly. 

One may now impose the condition q\q 2 = 1 and show that in non- vanishing cases 
the operator in (A. 4) is just a factor of J times that of (A.l). The first step is to 
substitute q 2 = 1/qi in Oi, and use the relative position index a = I — k to rewrite the 
double-summed expression in (A. 4) as 

j 

O l = J2( J - a )?2 TT((j)Z a ^Z J - a ). (A.5) 

a=0 




Figure 13: A circular array to indicate the J 
Z-fields and the impurities <j), tfj at positions k, I 
from a chosen origin. 



- 35 - 



We have used the fact that, for fixed a, there are J — a values of the original index k 
which make identical contributions to 0\. The sum over a in (A. 5) actually stops at 
a = J — 1, but it is useful to add the vanishing entry as we will see. The operator 2 is 
handled similarly using b — k' — I' as the relative phase index. By symmetry one finds 

j 

2 = ^(J-%?Tr(^V J - 6 ) 

6=0 
J 

= 5>g 2 a Tr(0Z>Z J - a ), (A.6) 

a=0 

where we have redefined a = J — b, used cyclicity and qi = l/q 2 in the last step. We 
see that the sum of (A. 5) and (A.6) is equal to J times the original on-shell operator 
(A.l) as claimed. 

We have gone into considerable detail in the simple case of 2 impurities in order to 
avoid an impossibly awkward discussion in the general case which we now outline. In 
a set of I > 3 impurities, repetition of the fields 0, ip occurs. However these fields are 
effectively distinguished in the construction of the operators because they carry different 
phases qj. Wick contractions in the correlators will then impose Bose symmetry. 

We therefore conceive of a set of I independent impurity fields i? % = 1, 
placed at arbitrary interstitial sites in the circle of figure 13, with assigned phase q k \ 
The analytic expression for the corresponding "off-level-matching-shell" operator is the 
sum I\ terms, one for each permutation of the impurities. The non-vanishing on-shell 
operator contains (/ — 1)! terms including all non-cyclic permutations. 

For 1 = 3, the first of six terms can be written as 

01= q\q l 2 q^^{Z k ^ x Z l ~ k ^ 2 Z m ~ l ^Z J - m ). (A.7) 

0<k<l<m<J 

There are two similar terms, 02,03 for the cyclic permutations (02,03 ? 0i) and 
(03,01,02) of the impurity fields, and three more for anticyclic permutations. With 
due diligence one may repeat the argument above for the case 1 = 2 and show that 
the sum Oi + 2 + O3 vanishes unless the level-matching condition qiq 2 q3 = 1 holds. 
The same property holds separately for the sum of the three operators for anti-cyclic 
permutations. The on-shell operator is the sum = O c + O a with the cyclic term 

a+b<J 

Oc= q a 2 q a 3 +b M<PiZ a <p2Z b <P 3 Z J - a - b ) (A.8) 

0<a, 0<6 
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and an analogous expression O a for the anti-cyclic permutation (0i,0 3 ,0 2 ), namely 

a+b<J 

o a = q a 3 q a 2 +b M^z a hzV z - a - b )- (A.9) 

0<a, 0<6 

With care, and with the 1 = 2 case as a model, one can insert the on-shell condition 
qi = l/(<72<?3) in the operator Oi, introduce relative position indices a — l — k,b — m — l, 
and rewrite 0\ as 

a+b<J 

1= Y (J-a-b^q^TrifaZOfcZOfoZ^-O). (A.IO) 

0<a, 0<6 

Cyclic symmetry implies similar expressions for 2 , 3 , namely 

c+d<J 

02= (J-c-d)qlq c 1 +d Tr( ( t )2 Z c (f )3 Z d ( t) 1 Z J - c - d y, 

0<c, 0<d 
e+f<J 

°3= (J-e-fK^MfoZ^ZffoZJ-^f). (A.ll) 

0<e, 0</ 

With the redefinitions c = b,d = J—a—b, f — a,e — J—b—a, creative use of the relation 
?i?2?3 = 1, and cyclic symmetry, one may show that the sum 0\ + 2 + O3 = JO c . 
The analogous result relating the on-shell sum of three anti-cyclic permutations to J 
copies of the anti-cyclic O a may be derived in the same way. This discussion shows 
that when the level matching condition holds, the off-shell operator for 3 impurities is 
just J times the on-shell operator O = O a + Ob- 

The 2-point correlation function of the operator with 1 = 3 impurities may be 
denoted by {(O c (x) + O a (x))(O c (0) + O a (0))). Planar diagrams come only from the 
diagonal terms (O c (x)O c (0)) and (O a (x)O a (0)}. Contributions from both terms are 
required for planar orthogonality in the two independent momenta. 

To order 0{g\ M ) we obtain 

j j 

(O(x)O(0)) = F(x) [l + g 2 YM NP(q z ,r t )ln\x\] J> 2 r 2 )' 5> 3 r 3 )"\ (A.12) 

Z=0 m=0 

P(q h n) = (1 - gf 1 )(l - ff 1 ) + (1 - g 2 )(l - f 2 ) + (1 - g 3 )(l - f 3 ). 

F{x) = (N/47i 2 x 2 ) J+3 is the product of J + 3 scalar propagators and the corresponding 
color factor. P{q i ,r i ) is of order 1/J 2 exactly like the 0(gy M ) expression for the two- 
point function of operators with two insertions (A. 2). In the expression for P{qi, r«) we 
have neglected higher order terms in 1/J that multiply terms of order 1/J 2 . 
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Figure 14: Cancellations at order g\ M . 



It is easy to understand how the structure of the interaction term arises. We are 
interested in planar contributions so we consider the nearest neighbor interactions only. 
The interactions can take place between any of the marked fields <pi and the neighboring 
fields Z. The interaction in the conjugate operator O has to take place between the 
same marked field and neighboring fields 2. Thus, we get the sum of three contributions 
with different phase dependence for each of the marked fields. 

We hope that the discussion for 7 = 2,3 impurities in this section makes the 
construction of the case of arbitrary / clear. 

B. Irrelevance of D-terms 

In this appendix we will show that the F-term interactions studied in the main text are 
the only interactions which need to be considered at order g\ M in Yang-Mills pertur- 
bation theory. The sum of the other contributions, from D-terms, gluon exchange and 
self-energy insertions, precisely vanishes at this order in both 2- and 3-point functions 
(for an example see figures 14 and 15). This simple but useful fact can be proved by 
minor modifications of the techniques used for the same purpose in the studies of BPS 
operators in either [7] or [54]. We use the technique of [7]. 
We are concerned with the operators 



and we will show that the sum all non-F-term contributions vanishes term-by-term in 
the expansion in phases of (OjOj) and {OjOj 1 Ty(2 j ' 2 )). 

The first relevant observation comes from inspection of the form of the D-term 
potential in the M = 4 Lagrangian which is 



V D = g 2 YM Tr ([Z, 2] [Z, 2} + 2[0, 0] [Z, 2} + 2ty>, $ [Z, 2} + 2[0, 0] ty, , (B.2) 



j 




(B.l) 
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Figure 15: Gluon exchange contributions to a three-point function. 



where we have dropped similar 4 and ip 4 terms which do not contribute to the cor- 
relators listed above. One now sees that all quartic vertices contribute to Feynman 
diagrams with the same combinatorial weight, independent of SU(3) flavor. Similar 
remarks hold for gluon exchange diagrams. The 1-loop self-energy insertion is also 
flavor blind. Thus for the purposes of this Appendix, each summand in (A.l) can be 
replaced by Tr(Z J+2 ). One can now simply use the result of [7] which shows that all 
order g\ M radiative corrections to (Tt(Z j+2 )Z j+2 ) cancel. Nevertheless, we will repeat 
the argument of [7] briefly because we will make a somewhat new application to 3-point 
functions below. 

The first step is to observe that gluon exchange diagrams and those from Vd ~ 
Tr([Z, Z] 2 ) have the same color structure, and must be summed over all pairs of lines 
in the second Feynman diagram of figure 14. Self-energy insertions on each line must 
also be summed. The following identity holds for any set of matrices M i? N: 

n 

Tr (M x • • • [M u N] ■ ■ ■ M n ) = 0. (B.3) 

i=l 

Let k = J + 2. Each diagram in figure 14 includes a sum over k\ permutations of the 
fields in Tr(Z fc ) relative to a fixed permutation of the fields of Tr(Z fe ). Let T ai ■ ■ ■ T ak de- 
note the fixed permutation of color generators of the fields in Tr(Z fc ), and let T bl ■ ■ ■ T bk 
denote one of the permutations of fields in Tr(Z h ). For each pair of fields % ^ j the 
gluon exchange or D-term has a color structure which may be expressed as a commu- 
tator with the generators T % and T J in the product T bl ■ ■ ■ T bk . Summing over all pairs, 
we obtain the net contribution 

k 

-fi(x)Tr(T ai ---T afc ) Tr(T h ■■■[T bi ,T e ]---[T b \T e ]---T bk ), (B.4) 
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Figure 16: Two representations of a diagram which cancels with others. 

where B(x) = b + bi \n(x 2 M 2 ) is the space-time factor associated with the interaction. 
We now use (B.3) on one of the commutators to rewrite (B.4) as 

k 

B(x) Tr(T ai ■ • ■ T ak ) ^ Tr(T bl ■ ■ ■ [[T\ T e ], T e ] • ■ ■ T bk ) 

i=l 

k 

= NB(x) Tr(T ai ■ • ■ T ak ) ^ Tr ( T&1 ■ ■ -T h ■ ■ -T bk ). (B.5) 

i=l 

In the last step we recognize [[ ,T e ],T e ] as the SU(N) Casimir operator in the 
adjoint representation which gives [[T a ,T e ],T e ] = NT a for any generator. The final 
sum thus has k identical terms. Each self energy insertion also contains the adjoint 
Casimir and has the form NA(x) = N(ao + a\ ln(x 2 M 2 )), and there are k such terms. 
The sum of all diagrams in figure 14 is thus 

kN(B(x) +A(x))Tr(T ai . . . T afc ) Tr(T fcl ■■■T bk ) ) (B.6) 

which must be summed over all permutations {5j = cr(aj) : % — 1 • • • k} and finally 
contracted on pairwise identical color indices. All manipulations above are valid for the 
case k = 2 which is known to satisfy a non-renormalization theorem. Hence B + A = 0, 
and radiative corrections (other than from F-terms) cancel for all k. Figure 16 shows 
a D-term diagram which cancels with others despite the intuition that a gauge theory 
vertex at the string interaction point, i.e. the saddle point of the toroidal stringy 
diagram, should be significant. 

Next we study the 3-point function {OjOj 1 Tr(Z j2 )) with J = Ji + J 2 . The flavor 
blind property again means that the JJ\ summands are all identical. Gluon exchange 
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interactions among pairs of lines are indicated in figure 15. Quartic vertices from Vd are 
similarly summed. We consider a fixed permutation of fields in Oj x and in Tr(Z j2 ) and 
sum over permutations (i.e. orderings of generators T Cl • • • T CJ + 2 ) in the central operator 
Oj and sum over permutations of fields from Oj 1 and from Tr(Z j2 ). For interacting 
pairs which are connected to Oj 1 the previous argument applies mutatis mutandis. 
Radiative corrections cancel when self-energy insertions on J\ lines are included. Idem 
for interacting pairs connected to Tr(Z j2 ). The remaining pair interactions include 
one line connected to each operator. For these we use the color structure to place one 
commutator in each position in Tr(T ai • • • T aj i ) and one commutator in each position 
within T bl ■ ■ ■ T bj z . The resulting structure is then 
Ji,J 2 

C(x, y, z) Tr ( Tai ' ' ' [ T *> T 1 • • • T<1Jl ) Tr ( Tbl " " " t Ti ' T 1 • • • Tbj2 ) Tr ( TC1 ' ' ' T CJ+2 ), 
i=i,j=i 

(B.7) 

where C(x,y,z) is a spacetime factor which need not be specified. However, for each 
fixed i, the sum on j vanishes by (B.3), and our task is complete. 

C. Feynman diagrams and combinatorics 

In this appendix we give a self-contained approach to the two-point correlation func- 
tion discussed in the text. The purpose is to provide the detailed basis of results for 
planar and genus one contributions and to outline an algorithm to calculate genus 
h > 1 results. We hope that the treatment below is readable both by physicists and 
mathematicians. 

Let us summarize the results of this appendix. We show that genus h two-point 
function in the free case is given by a sum of (Ah — l)\\/(2h+ 1) terms that correspond 
to the types of genus h Feynman diagrams with Ah nonempty groups of edges. The two- 
point function with a single interaction equals (27r) 2 (nn' + (n — n') 2 ) x (free case) plus 
sum of (Ah — l)(Ah — l)!!/[3(2/i + 1)] terms given by explicit formulas. The expression 
(27r) 2 nn'x (free case) comes from nearest neighbor interactions, the expression (2n) 2 (n— 
n') 2 x (free case) comes from semi-nearest neighbor interactions, and the remaining 
expression comes from non-nearest neighbor interactions. The latter correspond to the 
types of genus h diagrams with Ah — 1 nonempty groups of edges. For genus h — 1, 
there is exactly one type of diagrams with 4 groups of edges, and there is one type of 
diagrams with 3 groups of edges. For genus h = 2, there are 21 type of diagrams with 
8 blocks and 49 types of diagrams with 7 blocks, etc. 

It is natural to assume that, for 2 interactions, we need to go one level lower, i.e., 
the two-point function should be given by terms that depend on the free case and the 
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single interaction case plus new additional part given by diagrams with Ah — 2 groups 
of edges, and so on for higher h. This gives a natural hierarchy of genus h Feynman 
diagrams according to their numbers of blocks. 

It is important to note that the notation used in this Appendix differs from that 
used in the main body of the paper in some respects. In subsections 3.3 and 4.2 of 
the main paper we have used the symbols U = l-nm and V = 1-nn to parameterize the 
phases e tU and tV in BMN operators. In this appendix, we will usually use the symbols 
u, v for these quantities. The indices n, m will be denoted n, n'. The symbol A n>n r 
used repeatedly in sections 3 and 4 of the main body of the paper is identical to the 
symbol n') in this appendix. Finally, the contribution of non-nearest interactions 
denoted B nm in (4.10) is identical to G\ n {n, n') in the appendix. 

C.l Correlation functions 

For a positive integer J and two integers n and n 1 ', the operators J n and 3 n , are given 
by 

0£ = X>TWZtyZ J -') and 6 J n , = £ f v Tr(Z p jZ J - l 'fi , (C.l) 

1=0 l'=0 

where q = exp(27rin/J) and r = exp(27rm'/ J). (Here % = v^T). All fields Z, Z, 0, 
0, ip, ip are given by N x N Hermitian matrices. We will discuss the free two-point 
correlation function: 

{Ot(x)dl,{y)) = (4n 2 (x-y) 2 )- J - 2 Aj, N (n,n>) (C.2) 

and the two-point function with one interaction: 

(Ofa) [Z^\ [Z^}Oi(y)) = - 92Nl ^ X 2 ' yl Gj, N (n,n'). (C.3) 
The functions Aj^ N and Gj^ can be written as series in powers of N: 

A JtN = N J+2 J2 N ~ 2hA J and G J>N = N J+3 J2 N ~ 2hGh J- ( C - 4 ) 

h>0 h>0 

This is called the genus expansion of correlation functions, because A 1 } and G 1 } are given 
by sums over Feynman diagrams of some type drawn on an oriented genus h surface. 
As we will see, A 1 } is of order J 4/l+1 and G 1 } is of order J 4/l_1 as J — > oo. Let 

A h (n,n') = lim A h j/J 4h+1 and G h (n, ri) = lim G^/J^' 1 . (C.5) 

J— »oo J— >oo 
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The limits A(g 2 ;n,n') = \imAj, N /(J-N J+2 ) and B(g 2 ;n,n') = lirnGj^/U" 1 ■ N J+S ), 
where J, N — > oo so that J 2 ~ g 2 N for a fixed constant g 2 , can be written as 

A(g 2 ; n, n') = $> 2 ) 2/l A h (n, n') and G(^ 2 ; n, n') = $> 2 ) 2A G h (n, n>). (C.6) 

h>0 h>0 

For any integer values of n and n', A(g 2 ;n,n') and B(g 2 ;n,n') are analytic functions 
of #2- In particular, 

(fe ' )_ ^ ' (C.7) 
A(g 2 ; n, 0) = A(g 2 ; 0, n') = G(^ 2 ; 0, 0) = G(g 2 ; n, 0) = G(g 2 ; 0, n') = 0, 

for n, n' ^ 0. We will see that 

A h (n,n') = sum of {Ah + 1) • integrals, 

G h (n, n') = (2tt) 2 • {-nn' + (n - n') 2 ) • A h {n, n') + G^ n (n, n'), (C.8) 

where G^ n {n,n') = sum of 4 • (4/i - 1) • ^"^j" integrals. 

Each of these integrals is given by an explicit formula. As an example, we will present 
closed formulas for A h (n,n') and G h (n,n') for small values of genus h. In general, the 
expressions for A h (n,n') and G h (n,n') have the form 

polynomial in n and n' , . 

(C.9J 



(n ■ n') a ■ (n — n') b ■ (n + n') c 
C.2 Free two-point function via permutations 

Feynman diagrams for the free two-point function are basically given by permutations. 

Let us recall a few basic facts about permutations. A permutation of order m is 
a bijective map w : {1, . . . , m} — > {1, . . . , to}. Multiplication of permutations is given 
by composition of maps. All permutations of order to form the symmetric group S m . 
A cycle in a permutation w is a subset of the form {w(i),w 2 (i), . . . ,w r (i) = i}. In 
particular, each fixed point w(i) = i is a cycle of size 1. Thus each permutation gives 
a decomposition of {1, . . . , to} into a disjoint union of cycles. The number of cycles of 
w is the total number of cycles in this decomposition. Let the long cycle c G S m be the 
permutation that consists of a single m-cycle given by c : i i— > i + 1 (mod to). 

Feynman diagrams that describe Wick couplings in the free case (0^(x) 0^,(y)) 
are given by permutations w of order m = J + 2 corresponding to mappings between 
the fields in O^(x) and the fields in 0^,(y). The diagram with permutation w produces 
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a term with some power 

NC( W ) of the rank N of the 

gauge group U(N). The exponent 
C(w) is related to the genus h(w) of the corresponding diagram by Euler's formula: 

C(w)=m-2h(w). (CIO) 

In physics language, C{w) is the number of closed quark loops in the ribbonized diagram 
(a.k.a. fat graph) corresponding to the Feynman diagram with U(N) adjoint fields. 
Combinatorially, C(w) is the number of cycles in the twisted permutation c~ x w~ x cw G 
S m , where c is the long cycle in S m . The expression Aj^, which gives the free two-point 
function, can be written as the following polynomial in the variables q, f, and N 

m 

a j>n = Yl nC{w) E^" 2f " w " 2 ' ( c - n ) 

weS m ,w(l)=l i>2 

where the sum is over all permutations w such that w(l) = 1 (the first marked field 
is always contracted with 0), and the product is over alii ^ 1 (choice of position of 
the second marked field ip) . 

Let us say that the number h(w) — (m — C{w))/2 is the genus of a permutation 
w G S m . It is always a nonnegative integer because the parity of the number C(w) of 
cycles in the twisted permutation c~ 1 w~ 1 cw is the same as the parity of m. Clearly, 
the genus of any cyclic shift wc r of w is the same as the genus of w. Thus without loss 
of generality we will assume that w(l) = 1. The h-th term A 1 } in the genus expansion 
of A J>N is given by the terms in (C.ll) with w of given genus h. Since the only genus 
permutations are the identity permutation and its cyclic shifts, we have 

A°j = (qf)° + (qff + • • • + (qf) J . (C.12) 

In the next section we will show that A 1 } is a polynomial in q and f given by a sum of 
order J Ah+l terms q l fK 

C.3 Block-reduction of permutations 

For large values of J, the expression (C.ll) involves a summation over permutations 
of large orders. This makes it difficult to calculate the limit J — > oo of this expression. 
Nevertheless it is possible reduce a permutation w of arbitrary large order m and small 
genus h to a permutation a of order < Ah. The permutation o is the block-reduction 
of w and is formally defined below. 

For a permutation w G S rn with w(l) = 1, let us subdivide the set {l,...,m} 
into blocks as follows. If there are consecutive indices i, % + 1, . . . , % + k in the cyclic 
order such that w(i) = j,w(i + 1) = j + 1, . . . , w(i + k) = j + k then we combine 
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the indices i, i + 1, . . . , i + k into a single block. Here all indices i and values w(i) 
are understood modulo m. We will choose the blocks to be as maximal as possible. 
Thus the permutation w gives a subdivision of the set {1, . . . , m} into disjoint union 
of blocks Bi, . . . , Bk (where 1 £ formed by cyclically consecutive elements and a 
permutation a E of blocks. We call the permutation a the block-reduction of w. A 
permutation cr £ S^, > 2, is a block-reduction for some u> if and only if 

<t(1) = 1; <t(z + 1) ^ o-(z) + 1, for i = l,..., A; -1; and a{k)^k. (C.13) 

Let us say that a permutation cr 6 ^ is block-reduced if it satisfies these conditions; 
and let BFi£ be the set of block-reduced permutations of order k and genus h. We will 
assume that BRq contains one element — the empty permutation of order — which is 
the block-reduction of the identity permutation of any order. 

Figure 17 shows two genus h = 1 Feynman diagrams given by permutations 
1 2 5 6 7 3 4 8 9 10 and 1 2 9 10 6 7 8 3 4 5 11 12. The first permutation has 3 blocks B x = 
{8,9,10,1,2}, B 2 = {3,4,5}, B 3 = {6,7} and its block-reduction is 13 2. The 
second permutation has 4 blocks Bi = {11,12,1,2}, B 2 = {3,4}, B 3 = {5,6,7}, 
-B4 = {8, 9, 10} and its block-reduction is 14 3 2. 

Lemma C.l A permutation w and its block-reduction o have exactly the same genera 
h(w) = h(o~). For each block-reduced permutation a e Sf. of genus h, we have k < 4h. 

Proof - Recall that genus of w was defined in terms of the twisted permutation 
w = c~ 1 w~ 1 cw. The similar twisted permutation a for the block-reduction a of w is 
obtained from w by removing all its fixed points. Thus a and w have same genera. 
Also, a G Sk is a fixed-point free permutation. Thus its number of cycles is C(a) < k/2 
and the genus of a is h(cr) = (k - C(o))/2 > k/A. □ 

There is only one genus block-reduced permutation and two genus 1 block- 
reduced permutations 132 and 1432. For genus 2, there are 21 elements in BRg and 49 
elements in BR^. In general, the following statement holds. 
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Proposition C.2 For h > 1, the numbers of elements in BR 4/l and in BR 4/l _ 1 are 
given by* 

|BR 4 \| = and IBR.VxH^^-IBR^I. (C.14) 

Proof - Elements of the set BR 4h are in one-to-one correspondence with gluing of 
the 4/i-gon into a genus h surface. An element a G BR 4/l is determined by its twist 
a = c^ x o v ca G Sm, where c is the long cycle in S^. The permutation a is a fixed- 
point free permutation with 2h cycles. Thus o is given by the product of 2h commuting 
transpositions. Let us label the sides of the 4/i-gon by the numbers 1, . . . ,Ah (in the 
clockwise order) and glue the side labelled i with the side labelled cr(i), for % — 1, . . . , Ah. 
This will produce a genus h surface. 

The numbers of gluings of any 2/c-gon into a genus h surface were calculated by 
Harer and Zagier [51]. In the case k = 2h, their result implies the number of gluings is 
(Ah- l)!!/(2/i+l). 

For a' G BR 4 l ^_ 1 , the twist a' G S^-i should be a permutation given by the product 
of 2h — 2 commuting transpositions and a single 3-cycle. Such a' can be obtained from 
a, for a G BR^, by merging the first vertex 1 with the last vertex Ah and replacing 2 
transpositions (1,<t(1)) and (Ah,a(Ah)) with a single 3-cycle (I, a (Ah), a (I)). We will 
get all a' such that the vertex 1 belongs to the 3-cycle. In order to get all possible a' 
we need to take all cyclic shifts, which gives the factor Ah — 1, and then divide by 3, 
because we counted all elements 3 times. Thus IBR^J = (Ah — l)/3 • jBR^J. □ 

In order to recover a permutation w G S m from its block-reduction a G 5^ we need 
to know the sizes of the blocks Bi and the position of 1 in the first block, which is 
the placement of the marked field <fi in (C.l). This information can be encoded as the 
sequence (pi, ... , bk+i) of integers, where bi = \Bi\ for % — 2, . . . , k; bi is the number of 
elements of B\ after 1 and bk+i is the number of elements of B\ before 1 (in the cyclic 
order). So b\ + bk+i + 1 = \Bi\. This sequence satisfies the following conditions: 

&! + ••• + h+i = m - 1 = J + 1; h, b k+1 > 0; b 2 , . . . , b k > 0. (C.15) 

For a block-reduced permutation a G S k and a sequence bi, . . . , b k +i as above, let 
b[,b' 2 , . . . , b' k+1 be the sequence given by b\ = b a -i (i) fori = 1, . . . , k and b' k+1 = b k+1 . 
Then the h-th term in the genus expansion of the free two-point function can be written 

*Rccall that {Ah - 1)!! = 1 • 3 • 5 • • • {Ah - 1). 
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as 

fe+i 

•45 = Z E E 

k>l o-pRR' 1 i=l 

" fe v (C.16) 

^ ^+-+^-if 6 ' 1+ -" +6 -«-i(( ? f)° + (gf) 1 + • • • + (grf ') J , 

v 6i+-+6 fc+ i=J+l / 

where the sum is over all block-reduced permutations a of arbitrary orders and fixed 
genus h and the internal sum in the parenthesis is over sequences b±, . . . ,bk+i that 
satisfy conditions (C.15). 

The number of terms g*f J (each of absolute value 1) in the internal sum is of order 
jfc+i Thus only the terms with maximal possible value k = 4h survive in the limit 
J — > oo and the whole expression is of order J Ah+l . Let 

Ah 

A\n,n')=\m Ll ^- 1 . (C.17) 

Then A h (n, n') is given by the sum over a G BR^ and % = 1, . . . , Ah + 1. In the limit 
J — > oo, each sum over bi, . . . , bu+\ turns into a (k + l)-dimensional integral. In the 
next section we show how to compute integrals of this type. 

C.4 Calculation of integrals 

Let us fix r numbers u±, . . . , u r & C. For a positive integer J, let = e Mj//J , j = 1, . . . , r, 
and let 

Sj( qi ,...,q r )= ( C - 18 ) 
bi+---+6 r =j+i 

where 6i, . . . , 6 r > run over all decomposition of J + 1 into a sum of positive integers. 
We are interested in the asymptotics of Sj as J — > oo. Clearly, Sj is a sum of order of 
J r_1 monomials in the qj. The expression Sj/ J r ~ x is just the J-th Riemann sum for 
an (r — l)-dimensional integral. Thus 

F r ( Ul ,...,u r ) = lim S J {q u ...,q r )/J r - 1 = [ e UlXl+ - +UrXr d Xl ■ ■ ■ dx r . u (C.19) 

there the integration is over the (r — l)-dimensional simplex 

A r _! = {(xi, . . . , x r ) | Xi + • • • + x r — 1, x±, . . . , x r > 0}. (C.20) 
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Proposition C.3 The functions F r (u\, . . . ,u r ) are recursively determined by the fol- 
lowing relations. If U\ = ■ • • = u r = u then 

F r (u,u,...,u) =e u /(r- 1)!. (C.21) 

If Ui 7^ Uj for some i and j then F r is obtained though F r ^ x by the divided difference 
operator: 

F r -i(. . . , Uj, . . . , Uj, . . . ) — F r -i(. . . , Uj, . . . , Uj, . . . ) 

r r [U\, . . . , u r ) = , [y,.ZZ) 

Ui — Uj 

where Ui means that the variable Uj is omitted. 



Proof - In the first case u\ — ■ ■ ■ — u r — u we just integrate the constant e u over 
the simplex A r _i, whose volume is l/(r — 1)!. In the second case we may assume that 
i — r — 1 and j — r. The integral can be written as 



(C.23) 



where a stands for 1 — x\ — • rc r _ 2 - The 1-dimensional integral in the parenthesis is 

equal to (w r -i ~ u r )~ 1 (e UrCt — e Ur - ia ), which gives the right-hand side of the recurrence 
relation. □ 

For nonnegative integers a±, . . . ,a r , let 

F( ai ,...,a r )(ui, ...,u r ) = F fc («i, . . . ,u u u 2 ,. . . ,u 2 ,. . . ,u r ,. . .,u r ), (C.24) 

where k = a\ + • • • + a r and we have a« copies of U{, i — 1, . . . , r, in the right-hand side. 
The function F^ ai; ^ >ar ^(ui, . . . ,u r ) is invariant under simultaneous permutation of a^'s 
and UiS. Proposition C.3 gives the following recurrence relations for these functions: 

p p 

17 (...,di,...,aj— 1,...) (...,ai—l,...,a,j,...) _ _ \ i 

-T(ai,...,a r ) - It Ui T Uj, ai,aj > 1, 

F( ai ,...,a r ){ui, ...,U r )= F(... ia . +0j . v .. )aj ,...)(. .. ) if Mi = Uj, 

F(a u ...,a r ){ui, ■ ■ ■ ,U r ) = F(... ia . v . .)(..., if Oi = 0, 

F (a) («)=eV(a-l)!. 

where and -Uj means that the corresponding terms are omitted. 

The next theorem presents an explicit expression for all these functions. 



(C.25) 
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Theorem C.4 The function F r (ui, . . . , u r ) is a continuous function of ui, . . . , u r de- 
fined on C r . If Ui ^ uj for all i and j then F r is given by 

r 

F r ( U1 , ... )Ur )=Y,^ - u j) _i - ( c - 26 ) 

i=l j^i 

For arbitrary ui, . . . ,u r , the function F r is obtained from this expression by continu- 
ity. Also, F r with repeated arguments can be obtained by differentiation of the above 
expression as follow. For distinct ui, . . . ,u r and a±, . . . , a r > 0, we have 

i ? (a 1 +l,...,a r+ l)(Ml,...,Mr) = ' ' F r { Ul , . . . , U r ) . (C.27) 

a\\ a r \ 



Proof - The first claim follows from Proposition C.3 by induction on r. In order to 
prove the second claim, remark that the (k — l)-dimensional integral for F( i+i,...,a r +i)) 
k = r + ai + ■ ■ ■ + a r , can be reduced to the following (r — l)-dimensional integral: 



F (oi+ i,... )ar+ i)(«i,...,u r ) = / ( n ~7 J 

J*r-i \ i=1 ai! J 



Oil J 

Now this integral for F(ai+i,...,or+i) is obtained from the integral for F r = F^...^) by 
applying the differential operator n[=i( a «0~ 1 ('9/ , 9' u j) ai to the integrand. □ 

Formula (C.26) says that F r (u\, . . . ,u r ) is the top coefficient of the Lagrange in- 
terpolation of the exponent f(x) = e x at points ui, . . . ,u r . The second claim (C.27) 
can be reformulated in terms of the generating function, as follows: 

^2 Z T ■ ■ ' Z 7 F ( ai +l,...,a r +l)(ui: ■ ■ ■ ,U r ) = F r (tti + Z 1: . . . ,U r + Z r ) . (C.29) 

ai,...,a r >0 

For example, according to Theorem C.4, we have, for distinct u^v^w^ 
F {1) (u) = F (2) (u) = e u , F (3) (u) = e u /2 , 

i ? (i,D(«,«) = ^f, ^(2,i)(«.«) = £s-g^. (C30) 

7-1 / \ P u p" 

-T fl 1 in - "; l>, IU) = 7 V7 V — 7 T7 V + 7 T? \ • 

1 1 i (u — v)(u — w) (u — v)(v — w) (u—w)(v—w) 
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C.5 Formula for free two-point function 

In this section we put everything together and give a formula for free two-point function. 

For a permutation a G Sk and 1 < i < k + 1, let us define the numbers lb((x), 
lrj(o"), rlj(a), and rr; (a) as 

= #0' I j < <, ^0') < ^(0), lr *(^) = #0' I j < <, ^0') > (c 31) 

rliW = #{j | j > i, a'(j) < a'(i)}, rr^a) = #{j | j > i, a'(j) > a'(i)}, 
where a' G Sk+i is the permutation obtained from a by adding a fixed point k + 1. 

Theorem C.5 TTie term o/ £/ie genus expansion of the free two-point function is 
given by 

A H (n,n')= ^2 F (M^ + ^(a),Th(a),vu(a) + l), (C.32) 

where the arguments of all F 's are (27r^/— l(n — n'), 27r^/— In, — 27r^/— In', 0). 

The functions F are explicitly given by Theorem C.4. This formula is valid for any 
complex values of n and n'. It involves rational expressions in n and n' and in the 
exponents of — In, — 27r^/— In', and 2ti\/— l(n — ra'). If n and ra' are integers then 
all exponents are equal to 1. 

C.6 Example: free case, genus = 1, 2 

Assume that w = 27rn and t> = 2nn'. 

In the case of genus h — 1, there exists only one permutation <r = 1432 in BR4. 
Theorem C.5 gives the following expression for A l (n,n')\ 

A^n, n) = F (1)0 ,o,5) + ^(2,0,2,2) + ^(2,1,1,2) + ^(2,2,0,2) + ^(5,0,0,1) , ( c - 33 ) 

all functions F are in the variables (i(u — v),i u, —i v, 0). 

The above formula is valid for arbitrary complex n and n' . Assume now that n 
and n' are integers. Then e lu = e lv = 1 and A 1 (n, n') does not involve any exponents. 
In the cases when some of the arguments i(u — v), iu, —iv, and coincide, we can use 
the reduction relations (C.25) to simplify the expression (C.33): 

A 1 (0,0) = 5F (6) , 

A\n, 0) = ^(0, n') = F (5il) + F (4)2 ) + F (3 , 3 ) + F m + F (1 , 5) = 0, 

(C.34) 

A l (n,n) = 2F( 0j o,6) + ^(0,2,4) + ^(2,0,4) + ^(1,1,4), 
A l (n, -n) = F (1A 5) + 3F (2i2;2) + ^(5,0,1) , 
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where n is non-zero, the arguments of F's in the third line are (iu,—iu,0), and 
(2iu, iu, 0) in the fourth line. 

We can calculate the expression for A 1 (n, n') explicitly using Theorem C.4. There 
are several cases that depend on which pairs of arguments i(u — v), iu, —iv, and 
coincide. The function A 1 (n, n') is given by 



A\n,n') 



24 



1 I 1 I 7 
60 ~ 6v? ~ 

1 , 35 

12 m 2 8t? 



if n — n' — 0; 

if exactly one n or n' is 0; 

\in — n' are non-zero; 



if n 



-n are non-zero; 



; 1 , 2 \ -\ + 4 + 4 - — - , 2 otherwise. 



(C.35) 



The genus h = 2 free two-point function A 2 (n,n') can be written as a sum of 
9 ■ | BR 2 . | = 189 integrals. We can calculate all these integrals using Maple. Explicitly 
A 2 {n, n') is given by 



A 2 (n,n') 



21 
8! 



1 

2 4 5! 



43 
72-7! 



+ 



1 



504 u 2 

_ 1 
4-5! u 2 



+ 



10 u 



_ + , 143 =f 



12 m 6 



+ 



2 6 i 2J53 
3 2 m 4 3m 6 



if n — n' — 0; 

if exactly one n or n' is 0; 

n' are non-zero; 
if n — —n' are non-zero; 



(C.36) 



some polynomial of degree 18 in u and v otherwise 

u 4 v 4 (u— v)& (u+v) 4 

We have skipped the numerator in the last case. 

For genus h = 3, there are 11M/7 = 1485 elements in BR^ 2 and each gives 13 terms. 
In total we need to calculate 19305 terms. This can also be done on a computer. 

C.7 Two-point function with an interaction 

The Feynman diagrams that correspond to Wick couplings in the case of the two-point 
function with a single interaction (0(x) [Zip] [Zip] 0(y)) can also be easily described 
in terms of permutations. A coupling is given by a permutation w G S n of order 
m = J + 2 with w(l) = 1 together with a choice of two indices i,j — 2,...,m,i^j, 
that correspond to the interacting fields. The corresponding diagram is obtained from 
the free Feynman diagram, given by permutation w, by adding an "interaction edge" 
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i ~» j between the ith and jth edges of the free diagram. Each such diagram given 
by a triple (w,i,j) produces four terms in Gj^ that correspond to four terms in the 
expansion of the product of commutators [Zip] [Zip]. All these terms come with a 
factor N c ( w ' hj \ where the exponent C(w,i,j) counts the number of closed loops in the 
ribbonized Feynman diagram. It is not hard to express this number combinatorially in 
terms of cycles in permutations: 

C(w,iJ) = #{cycles in w} + S w(i)t (C.37) 

where w = t w ^ t j-ic^w^cw and t w ^j-i is the transposition of w(i) and j — 1. Notice 
that the number C(w, depend on the order of % and j. The genus of the Feynman 
diagram associated with triple (w, is combinatorially determined by Euler's formula 

h(w,i,j) = (n + 1 - C(w,i,j))/2. (C.38) 

The two-point function with an interaction Gj^ can now be written as 

G J>N = Yl N J+3 ~ 2h ^ x 

weS m ,w(l)=l i,j=2,...,n;ift (C.39) 

x [q i f w< ^ > — q i f w ^ — qif w( -^ + q^f w<y ^ (qf)~ 2 . 
The h-th term of the genus expansion is given by the sum 

G) = {(l ifW{i) - <l ifW(j) - <l j r w{€) + q j f w{j) ) {qf)- 2 (C.40) 

over all triples (w,i,j) with fixed genus h(w,i,j) = h. 
We will see that G 1 } is of order J 4/l_1 . Let 

G h {n,n') = Urn G h j/J 4h -\ (C.41) 

J^oo 

We have h(w,i,j) > h(w), where h(w) is the genus of the free Feynman diagram 
as defined in Section C.2. According to Section C.3, the total number of permutation 
w G S m of genus h is of order J 4h . This implies that the total number of all triples 
(w,i,j) with h(w) < h(w,i,j) — h is of order J 4( - h ~^J 2 = J 4h ~ 2 and each of these 
triples give 4 terms of absolute value 1. Thus, in the limit of expression (C.39), the 
pairs such that h(w,i,j) > h(w) will not make any contribution to G h (n,n'). 

It is natural to subdivide all possible choices for the interaction that does not 
increase the genus h(w,i,j) = h(w) into the following there classes: nearest-neighbor 
interactions (j = i + 1 and w(j) = w(i) + l), semi-nearest neighbor interactions (exactly 
one of the conditions j = i + 1 or w(j) = w(i) + 1 holds), and non-nearest neighbor 
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interactions (j ^ i + 1 and w(j) ^ w(i) + 1). Let Gjj e (n, n'), G^(n,n'), and Gj n (n, n') 
be the contributions to G h (n,n') of these tree case, respectively. Thus 

G h (n, n') = G h ne (n, n>) + G h sn (n, n') + G h nn (n, n>) (C.42) 

Let us show how to calculate these three expressions. First, for a given permutation 
w G S m , we describe all pairs (i, j) such that h(w, = h(w). Suppose that B ± , . . . , B k 
are the blocks of a permutation w G S m . Let us connect all blocks B 1: . . . , B k by directed 
edges of two types " — >" and " — ■>" as follows: 

i?! < — i? 2 < — • ■ • < — B k < — Bi 

R R ♦ d + o (C43) 

An alternating chain of blocks of length / is a chain of the following type: 

B S1 — » 5 S2 — > B S3 — » fi S4 — > B SB -■»... (C.44) 

The length of a chain is its number of edges. An alternating cycle of blocks is a closed 
chain B S1 ---> ••• — > B Sr , = B Sl . 

Lemma C.6 Let w G S m and let i,j G {2, . . . ,m}, i ^ j. Then h(w,i,j) = h(w) if 
and only if one of the following two conditions is satisfied: 

1. j = i + 1 and i,j belongs to the same block B s . 

2. i the last element in some block B s ; j is the first element in some block B t ; and 
B s is connected with B t by an alternating chain of odd length: 

Bs * B t . (C.45) 



According to Section C.3, the total number of permutations w G S m of genus h 
with k blocks is of order J k . Since G 1 } is of order J 411 " 1 , only permutations with k = 4h 
or k — Ah — 1 blocks can give a nonzero contribution to G h (n,n f ). 

Let us first consider the case k = Ah. We have two options: 
I. Nearest neighbor interactions: Suppose that % and j are such as in case (1) of 
Lemma C.6. Then 

qifw(i) _ gifw(j) _ qjfw(i) _|_ qjfw{i) — — g)(l — f) q i r w ^. (C.46) 

The contribution of these terms to G 1 } is asymptotically equal to (1 — q)(l — f) J 4h+1 Aj. 
Note that it of order J 4/l_1 , because (1 — — f) is of order J~ 2 . Thus the contribution 
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of these terms to G h (n,n') is equal to 

G h ne {n,ri) = uvA h {n,n'), (C.47) 

where A h (n, n') is the h-th. term in the free two-point function given by Theorem C.5. 
As before, u = 2im and v = Inn! . 

II. Semi-nearest neighbor interactions: This is case (2) of Lemma C.6. For 
k = Ah all possible alternating cycles of blocks have length 4. Thus the only possible 
pairs % and j are the following: % is the last element of B s ; j is the first element of B t \ 
and (B s ---» B t or B s < — B t ). In this case we can recombine the terms as follows: 

(gV^W _ qifu>ti) _ qifw(i) _|_ gif^O')) — 

(C.48) 





(2 - q - f) > gV w + (2 - g" 1 - f" 1 ) > g J f 



where LAST(u>) is the set of last elements in blocks of w and FIRST(w) is the set of first 
elements in blocks of w. The right-hand side of this expression involves 8h terms, each 
come with a factor (2 — q — f) ~ — \/— T(u — v) J -1 or (2 — q~ l —f~ 1 ) = >/— T(u — u) J -1 . 
Since the number of elements w with a G BR^ is of order J 4h , the total contribution 
of terms of this type is of order J 4/ V -1 = J 4/l_1 . In the limit J — > oo, the sum over 
u> of given genus /i becomes a finite sum integrals given by Theorem C.4. Using the 
notation of Theorem C.5, we can write two terms in the right-hand side of (C.48) as 
the following sums over a and i such as in (C.32): 

(2-q-r) J2 1 JfWU) = - v ) E ^+1^^)' 

jeLAST(w) a,i (C 49) 

(2 - H ' - r- 1 ) ^ fW(3) = " «) E ^A^+D. 

jeFIRST(w) <r,i 

where the arguments of F's are (v^— T(u — f), \/— Tit, — \f— Tv, 0). According to (C.25), 
the sum of the two summands in the right-hand sides of (C.49) is 

(u - vf F(u.+i )lr . )r i. )IT . + i), (C.50) 

which is exactly (u — v ) 2 times the summand in (C.32). Thus the total contribution of 
semi-nearest interactions is equal to 

G h sn {n, ri) = {u- v) 2 A h (n, ri), (C.51) 

where A h (n, ri) is the free two-point function given by Theorem C.5. If n — ri then 
G*(n, n) = 0. 
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Let us now consider the case k = Ah — 1. The total number of permutations w 
with a E BR^_ 1 is of order J 4,l_1 . In this case we may ignore the nearest neighbor 
interactions because their contribution comes with a prefactor (1 — q)(l — f) ~ J~ 2 , 
which makes the total order J 4 ' 1-1 J~ 2 subdominant to J 4h_1 . In this case all alternating 
cycles of blocks have size 4, except a single cycle of size 6. The 4-cycles produce semi- 
nearest interactions that come with a prefactor (2 — q — f) or (2 — q" 1 — f~ l ) of order 
J -1 , which again makes their contribution negligible in the limit J — > oo. 

III. Non-nearest neighbor interactions: Suppose that k = Ah — 1 and w is a 

genus h permutation with Ah — 1 blocks. Then there exists a unique alternating 6-cycle 

B S1 — > B S2 — > B S3 — > £? S4 — >■ £? S8 ---> £? S6 — >■ _B Sl . (C.52) 

Let f r and l r be the first and the last elements, respectively, in the block B Sr . There are 
only 3 possible choices = (h, f±), (I3, fe), (h, $2) for the interaction edge, whose 
contribution survive in the limit J — > 00. Each of these 3 pairs produces 4 terms, which 
gives the following 12 terms: 

qhy;w(h) _|_ qf4fw(j 4 ) _ qh^w(f 4 ) _ qf 4 ^w(h) 

qhfw{h) _|_ qhfw(k) _ qhfw{k) _ qhfw{h) (C.53) 
qhfw(i 5 ) _|_ qf2fw(f 2 ) _ qkfw(f2) _ qhfw{h) _ 

Since the number of genus h permutations w with Ah — 1 blocks is of order J 4h_1 , the 
total contribution of these terms to G 1 } is again of order J 4ft - _1 . The sum of the 12 
terms in (C.53) over all genus h permutations w with Ah — 1 blocks can be written as 
the sum 

(q 1 - 1 ^^ - 2q i f Mt) + q i+1 f^ i+1 ^ , (C.54) 

w 

over all genus h permutations w with Ah blocks such that w contains a block B sing with 
a single element, |-B S in g | = 1, and i 7^ 1 is the position of the block -B S i ng in w. 

Indeed, there are the following 3 ways to transform w with Ah — 1 blocks to a 
permutation w with Ah blocks by inserting a block with a single element and preserving 
the genus. Let be the permutation obtained from w by inserting a new edge to its 
Feynman diagram between the blocks B\ and Bq on the top and the blocks B3 and B4 
on the bottom. Similarly, we construct the permutation by inserting a new block 
with a single element between B 3 and B 2 on the top and B 5 and B 6 on the bottom; and 
the permutation by inserting a new block between B 5 and B A on the top and Bi 
and B 2 on the bottom. One can easily check that the sum of the summands in (C.54) 
for 3 permutations w = w^\w^ 2 \w^ produces exactly the 12 terms in (C.53). 
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As before, in the limit J — > oo, the sum (C.54) reduces to a finite sum of integrals 
over block-reduced a G BR^. Using the notation of Theorem C.5, we can write the 
contribution Gj n (n,n') of non-nearest neighbor interactions as follows: 

4ft 

G h nn (n,n>) = -2 £ £ F (IWri , rli , rri) + 
*eBR% h i=2 



+ ^ I ^ ^(114+1,^,^,^-1)+ ^(lli+l.lri.rli-l.rri) J + (C.55) 

ff £BR 4 \ y=l,..., 4h-l i=l,... ,4/i-l 



+ ^ 1 5^ ^(lli-l.lri.rli.rri+l) + ^(ll <) lr i -l,pl<,rr i +l) 
ctGBR^ \i=3,...,4/i+l j=3,...,4/i+l 

where the arguments of all F's are (\/— T(u — v), y/^lu, — >/— If, 0). 
C.8 Example: interaction case, genus = 1 

Then the contribution of nearest and semi-nearest interactions is equal to 
G]Jn, ri) + Gl n {n, ri) = (uv + (u - v) 2 ) A\n, ri) = {u 2 + v 2 - uv) A\n, ri). (C.56) 



This reduces to the free genus 1 case given in Section C.6. 

Formula (C.55) gives the following expression for G nn (n,ri): 

GL(n,n') = -2(F (1)2i0) i) + F (lil)1>1) + F (1)0)2> i)) + 

+^(1,0,0,3) + ^(2,0,1,1) + ^(2,1,0,1)+ (C.57) 
+-^(1,0,1,2) + ^(1,1,0,2) + ^(3,0,0,1), 

where, as usual, the arguments of all F's are (i(u—v),iu, —iv, 0). Using relations (C.25), 
we can simplify this expression: 

GU0,0) = GUn,0) = GU0,n') = 0, 

Gl n (n, n) = 2(F( ,o,4) - ^(0,2,2) - ^(2,0,2) + ^(0,1,3) + ^(1,0,3) - ^(1,1,2)), (C.58) 
G i n (n, -n) = F(i >0)3 ) + ^(3,0,1) + 2 -^(i,i,2) + 2F(2,i,i) - 6F (li2 ,i). 

where n is non-zero, the arguments in the second line are (iu, —iu, 0), and the arguments 
in the third line are (2iu,iu, 0). 
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Finally, using Theorem C.4, we obtain 

G\n,ri) = {u 2 + v 2 - uv) A\n,n') + G)Jn,ri), (C.59) 

where 



G)Jn,ri) = < 



if n or n' is 0; 



\ + iin — n' are non-zero; 



15 
~2u 2 



if n — —n' are non-zero; 



(C.60) 



— + i 2 ,i otherwise. 



D. Effective operator approach to Wick contractions 

In this appendix we outline an approach to handling the color combinatorics of planar 
and toroidal Feynman diagrams which makes use of perhaps more familiar techniques 
based on Wick contractions. To demonstrate these methods we consider the planar 
and genus one contributions to the two-point function at 0{gy M ). This provides an 
independent check of many of the calculations presented elsewhere in this paper. 
Once again the BMN operators are, 

j J 
O = Q l Tr {<f>Z l tjjZ J - 1 ) and O = ^ f 1 ' Tr (z J - l '$Z l 'fi) (D.l) 

1=0 l'=0 

with q = exp(2mn/ J) and r = exp(2mn' / J) . 

As discussed in Appendix B it is only necessary to consider the F-term interactions 
in the Af = 4 SYM action. These can be written, 

-4^ M Tr([Z,0p,0]) = 4g 2 YM r b Z^ b f^Z'tp (D.2) 

where the f pab are the structure constants of SU(N), and we trivially extend them to 
U(N) by adding the N x N matrix T° = I/y/N to the standard set of N 2 - 1 SU(N) 
generators T a , a — 1, ... , N 2 — 1. The full basis of U(N) generators is then normalized 
byt 

(ab) = Tr T a T b = 5 ab . (D.3) 

Using equation (D.2) we define an effective operator, O eff , as the sum of Wick 
contractions of and each Z in (D.l) with the factor f pab Z a (j) b of the interaction. After 

TThis normalization differs from that common in physics by a factor of 2 so that the f pab differ by 
V2- 
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X 



trivial manipulation this produces 
j 

O cS = -iJ2 ql 

1=0 

1-1 j-i-i 

^2Tr([T a ,T p }Z m T a Z l - m - 1 4>Z J - 1 ) + Tr([T a , T p \Z l ipZ J - l - m - 1 T a Z m ^ 

\m=0 m=0 

(D.4) 

From now on we compactify the notation by replacing all explicit generators by 
their index values, i.e. T a — > a and replace the explicit trace of an arbitrary N x N 
matrix M by Tr(M) — > (M). The following 'splitting/joining' rules can then be used 
to evaluate traces and products thereof which involve summed repeated color indices: 

(MaM'a) = (M)(M') (ab) = 5 ab 
(Ma)(aM') = (MM') (a) = VN5 a0 (D.5) 

aa = NI ( ) = N 

These follow from results derived in ref . [55] . 

Finally it is useful to write O e s = 0\ s + 0"l s where 0\ s includes only the 'nearest 
neighbor interactions i.e., the m = terms in (D.4). The m > terms are contained in 
Og ff and represent interactions between fields which are non-nearest neighbors. These 
operators may be expressed as 

j-i 

Ol s = -iN(q - 1) £y( P Zty£ J -'- 1 ); (D.6) 

1=0 

j-i 

g 

Q 



2 cS = -i-^ 5>Z m )(Z J - m -V)(l + q- 1 ~q m - q- m - 1 ) 



m=l 
J-1 J 

" i E E - ?- m - 1 )(^ m )(p^- m -V^')- (D.7) 

m=l l=m+l 

In obtaining equations (D.6) and (D.7) we have applied the U(N) trace identities in 
equation (D.5) and explicitly summed several of the geometric progressions in q l by 
reversing the order of summations. 

With equations (D.6) and (D.7) at our disposal we may now discuss the order g\ M 
contributions to the correlator 

(O{x)0{y)) = (Ol s Ol s ) + (Ol s 2 cS ) + (OlA) + (0 2 cS 2 cS ). (D.8) 

From now on we omit reference to the space-time points x, y since our main concern 
is to capture the order N J+3 planar and order iV J+1 genus one contributions to this 
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correlator. This is handled here by performing the Wick contractions of the fields Z, Z 
and evaluating the ensuing color contractions and traces. Evaluating the various terms 
in equation (D.8) requires the following 'contraction identities', 

Ti(Z a Z a ) = N a+1 + * ^ N a ~ x + 0(N a - 3 ) 

Tr{Z a ) Ti{Z a ) = aN a + 0{N a - 2 ) (D.9) 
Tr(Z a ) Ti(Z a+h Z b ) = a(b + l)iV a+6 + 0{N a+h ' 2 ) 

Tr(Z a Z c Z b Z d ) = (min(a, b, c, d) + 1) N a+b+1 + C(iV a+6 - 1 ) with a + b = c + d 

which are derived by counting the number of ways one may perform the Wick contrac- 
tions within each trace structure while obtaining a maximal power of N. The remaining 
0(N J+1 ) contributions are the semi-nearest and no n- nearest diagrams of section 4 of 
the main text. In the effective operator formalism the semi- nearest diagrams are given 
by the second and third terms in (D.8), and computation shows that they vanish in the 
special case q = r we are considering. The non-nearest diagrams are given by the last 
term in (D.8). It is straightforward to evaluate and sum the relevant Wick contractions 
and obtain 

j-i j-i 

(OlA) = N\q - \){q-\)Y^n{Z l Z\Z J - l -'Z J -~ 1 ^) (D.10) 

1=0 I=o 

for which planar contributions are only possible for 1 = 1. It is interesting to note that 
up to the (q — l)(q — l) prefactor this is exactly the same expression which generates the 
all-genus polynomial at order gy M . Use of the identities in equation (D.9) and explicit 
evaluation of the resulting sums yields an expression which in the limit of large J is 
given by 

(OlA) = g\A nn (1 - n 2 A' ln(AV)) , 

where A nn = ( ^- - \ — 1 + 7 ) . (D.ll) 
\60 6 {2nn) 2 (2mip J 

Notice that the non-planar corrections are of order g\ as has been emphasized through- 
out this paper. In order to obtain the remaining 0{N J+1 ) contributions we need only 
consider the final term in equation (D.8) since in the special case we are considering, 
i.e., q = r, the other terms conspire to precisely vanish at leading order in a large J ex- 
pansion. Evaluation of the final term in equation (D.8) is a straightforward calculation 
which yields, 

gjXHA^ n 5 \ 
\°«0«)- 47r2 U + 2 ^ ] (D - 12) 
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Putting the results of equations (D.ll) and (D.12) the complete planar and torus 
contributions to the g\ M piece of the two-point function is given by, 

(1 + S 1A,„) (1 - n'X' ln(AV)) + (I + JL^ (D.13) 

which is precisely the same answer as obtained in the main text. Extending these 
techniques to arbitrary q, r is straightforward and in all cases is found to agree with 
the methods used elsewhere in this paper. 
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